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1.0  INTRODUCTION 

This  theoretical  and  experimental  analysis  of  the  constant-area, 
supersonic-supersonic  ejector  was  prompted  by  current  interest  in  the 
high-energy,  chemical  laser  and  the  unique  gas  dynamic  problems  associated 
with  that  device. 

In  a typical  high-energy,  chemical  laser,  Figs,  1.0-1  and  1.0-2, 
hydrogen  and  flourine  are  precombusted  to  form  free  flourine  atoms  which, 
with  the  other  gaseous  products  of  combustion,  accelerate  into  the  laser 
cavity  through  a series  of  supersonic  nozzles.  Additional  streams  of 
secondary  gases  enter  the  laser  cavity  through  a series  of  supersonic  noz- 
zles to  form  alternate  interleaved  streams  of  precombustor  and  secondary 
gases.  These  streams,  initially  at  high  Mach  numbers  (3  to  7)  and  low 
static  pressures  (5  to  200  Terr)  mix  and  react  to  establish  the  lasing 
zone  by  chemically  producing  population  inversions  of  selected  species. 
Accompanying  these  chemical  reactions,  a significant  quantity  of  heat  is 
released  into  the  laser  cavity  flow  which  tends,  qualitatively,  to  in- 
crease the  static  pressure,  to  decrease  the  stagnation  pressure,  and  to 
decrease  the  Mach  number  of  the  "mixed"  supersonic  flow.  At  the  laser 
cavity  exit,  the  hot  (T  - 1500  K) , corrosive,  supersonic  (1.5  < M < 3.5) 
"mixed"  flow  at  low  pressure  (10  < P < 60  Torr)  must  be  "pumped"  by  some 
form  of  diffuser-ejector  system  to  atmospheric  discharge  conditions 
(760  Torr)  in  order  to  start  and  sustain  the  lasing  process.  A pumping 
system  wi.ich  ,s  ideally  suited  to  this  application  would  have: 

1.  A potential  compression  ratio  in  the  range  80  to  8, 

2.  Simplicity  of  design  with  high  resistance  to  hot,  corrosive 


gases, 


2 


3.  The  capability  to  isolate  the  laser  cavity  flow  from  minor  per- 
turbations in  the  downstream  conditions  and/or  pumping  system, 

4.  The  flexibility  in  operating  point  and  performance  to  allow  for 
variations  in  laser  performance  and  flow  conditions, 

5.  The  capability  for  short-duration,  transient  start-up, 

6.  Minimum  pumping  resource  requirements,  and 

7.  Compact  design  with  portability  and  mobility  as  goals. 

Current  high-energy  chemical  laser  designs  have  incorporated  the 

constant-area,  supersonic  diffuser  coupled  with  the  constant-area, subsonic- 
supersonic  ejector,  Fig.  1.0-3,  as  a pumping  system  since  the  operation  of 
this  system  is  fairly  well  understood  while  providing  the  most  obvious,  if 
not  satisfactory,  solution  to  the  above  requirements.  Nevertheless,  the 
quest  for  an  improved  pumping  system  continues  with  various  modifications 
of  the  conventional  diffuser-ejector  as  candidates  [1]. 

One  such  candidate  for  a high-performance  chemical  laser  pumping 
system  is  the  supersonic-supersonic  ejector,  Fig.  1.0-4,  the  subject  of 
this  investigation.  In  this  system,  the  diffuser  is  eliminated  and  the 
supersonic  laser  cavity  stream  is  pumped  directly  by  the  ejector.  Thus, 
the  supersonic-supersonic  ejector  offers  a potential  alternative  in  that 
the  desirable  characteristics  of  the  conventional  diffuser-ejector  are  re- 
tained with  a reduction  of  size  and  a possible  increase  in  performance. 

1.1  REVIEW  OF  PREVIOUS  WORK 

The  ejector  has  been  in  use  for  many  years;  indeed,  the  literature  is 
filled  with  a virtual  multitude  of  ejector  related  papers*,  the  author 


*See  APPENDIX  7.0. 
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having  compiled  over  300  entries  dating  from  1892.  With  the  exception  of 
three  supersonic  wind  tunnel  studies,  none  of  these  papers,  to  the  author's 
knowledge,  have  addressed  the  problem  of  pumping  a supersonic  stream 
directly  by  an  ejector,  which  is  not  particularly  surprising  since  most 
commercial  applications  involve  the  pumping  of  subsonic  or  nearly  stagnate 
streams. 

The  three  wind  tunnel  studies  [2 . 3 , 4 j were  part  of  an  experimental  in- 
vestigation to  study  the  effects  of  auxiliary  air  injection  on  the  pressure 
recovery  of  variable  geometry,  supersonic  wind  tunnel  systems.  In  each 
case,  the  auxiliary  air  was  injected  through  a supersonic  nozzle  at  the 
downstream  end  of  the  supersonic  test  section,  thus  in  essence  forming  a 
supersonic-supersonic  ejector.  The  first  two  investigations  [2,3]  showed 
that  the  resultant  pressure  recovery  with  injection  was  not  as  good  as 
those  attained  with  variable-geometry  diffuser  tunnels;  however,  the  later 
study  by  Hasel  and  Sinclair  [4]  demonstrated  a significant  improvement  in 
total  system  pressure  recovery  with  auxiliary  injection. 

The  general  methods  of  analysis  employed  in  existent  subsonic- 
supersonic  ejector  models  apply  equally  well  to  the  supersonic-supersonic 
ejector  provided  their  application  is  consistent  with  the  physical  phenom- 
ena. These  methods  of  analysis  are: 

1.  The  one-dimensional  analysis, 

2.  The  method-of-characteristics , 

3.  The  method  of  integral  relations,  and 

4.  The  finite-difference  method. 

In  the  one-dimensional  model  of  the  subsonic-supersonic  ejector  as 
executed  by  Fabri,  et  al.  [5,6],  the  conservation  equations  were  applied 
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| to  a control  volume  contained  within  the  ejector  mixing  tube  assuming  uni- 

form velocity  and  pressure  distributions  of  the  primary  and  secondary 
streams  at  the  tube  entrance  and  a uniformly  mixed  stream  at  the  tube  exit. 
A second  control  volume  was  used  to  predict  the  operation  of  the  ejector 
when  the  secondary  stream  choked  within  the  mixing  tube.  The  primary  and 
secondary  streams  were  assumed  to  remain  distinct  and  to  be  isentropic  from 
their  point  of  confluence  to  the  secondary  choking  location.  Consequently, 
the  condition  that  the  static  pressures  be  equal  at  the  boundary  between 
the  primary  and  secondary  streams  is  not  satisfied. 

A modified  or  quasi-one-dimensional  model  has  also  been  applied  to 
ejectors  with  short  mixing  tubes  m which  case  the  secondary  and  primary 
streams  are  assumed  to  remain  distinct  from  entrance  to  exit  [7]. 

Addy  and  Chow  [8-11]  developed  a more  sophisticated  approach  for  the 
prediction  of  secondary  stream  choking  within  the  mixing  tube.  In  this 
model,  the  secondary  stream  was  treated  by  the  conventional  methods  of  one- 
dimensional gas  dynamics  while  the  primary  flow  field  was  obtained  from  the 
two-dimensional  method-of-characteristics  for  steady,  irrotational , super- 
sonic axi symmetric  flows.  This  method  allows  the  static  pressures  to  be 
matched  at  the  boundary  of  the  streams  since  the  primary  stream  may  have  a 
nonuniform  pressure  profile;  although,  the  pressure  must  be  uniform  across 
the  secondary  stream.  The  simultaneous  solution  of  the  two  flow  fields 
satisfies  the  choking  criteria,  i.e.,  a Mach  number  of  unity  at  the  minimum 
flow  area,  was  then  corrected  for  viscous  effects  by  superimposing  the 
mixing  layer  on  the  inviscid  boundary  between  the  primary  and  secondary 


streams . 
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Hewlett  and  Chow  [12]  added  more  detail  to  the  Addy  and  Chow  model  us- 
ing the  method  of  integral  relations  for  computation  of  the  secondary  stream 
but  retaining  the  method-of-characteristics  for  the  primary  stream;  the 
static  pressure  was  matched  at  the  boundary  of  the  streams.  The  choking 
criteria  for  the  secondary  stream  was  developed  based  on  a singularity  in 
the  integral  relations  describing  the  secondary  flow. 

Hill,  et  ai.  [13],  applied  the  method  of  integral  relations  to  both 
the  secondary  and  primary  streams  but  later  [14,15]  adopted  a finite- 
difference  model  which  does  not  attempt  to  separate  the  secondary  and 
primary  flows. 

Each  of  the  models  discussed  has  its  own  particular  advantages  and 
disadvantages. 

The  one-dimensional  model  is  well  suited  to  broad-band  parametric 
studies  of  ejector  operation  since: 

1.  It  is  computationally  simple,  each  operati  ig  point  being  deter- 
mined by  the  direct  solution  of  a set  of  algebraic  equations; 

2.  It  applies  equally  well,  at  least  in  theory,  to  any  constant-area 
ejector  configurations;  and 

5.  It  is  quite  reliable  as  long  as  the  assumption  of  one -dimensional 
or  quasi- one-dimensional  flow  is  satisfied. 

The  one-dimensional  model  has  certain  disadvantages  in  that: 

1 It  is  restricted  to  the  constant-area  and  constant-pressure 
ejectors  where  any  pressure-area  surface  forces  acting  in  the 
flow  direction  are  eliminated  from  the  momentum  equation; 

2.  It  may  seriously  error  in  flow  regimes  which  are  highly  two- 


dimensional in  nature; 
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3.  It  provides  no  insight  into  the  actual  flow  phenomena;  and 

4.  It  requires  some  a priori  knowledge  of  ejector  operation  to  pre- 
dict any  limiting  conditions  such  as  choking  of  the  secondary 
stream. 

The  method-of-characterist ics , method  of  integral  relations,  and 
finite-difference  models  have  3 A the  advantages  of  sophistication  in  that: 

1.  They  may  be  applied  to  all  ejector  configurations  including 
variable-area  geometries; 

2.  They  provide  field  descriptions  of  increasing  detail; 

3.  They  satisfy  the  physical  condition  of  continuity  of  static  pres- 
sure across  the  boundary  between  the  primary  and  secondary  streams 
and 

4.  They  produce  good  results  over  all  phases  of  ejector  operation. 

On  the  other  hand , these  models  require: 

1.  Considerable  knowledge,  even  empirical  relations,  taken  from  prior 
experimentation  for  their  development;  and 

2.  Significant  amounts  of  computer  time  for  program  development  and 
convergence  problems  which  restricts  their  use  for  parametric 
studies . 

1.2  STATEMENT  OF  THE  PROBLEM 

This  theoretical  and  analytical  analysis  of  the  constant-area, 
supersonic-supersonic  ejector  was  conducted  to: 

1.  Develop  a simplified  mathematical  model  for  predicting  the 
operating  characteristics  of  a constant-area,  supersonic- 
supersonic  ejector  which  is  suitable  for  parametric  evaluations 
and  optimization  procedures; 
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2.  Provide  quantitative  and  experimental  data  for  verification  of 
the  theoretical  model  and  identification  of  problem  areas  not 
indicated  by  the  theoretical  analysis;  and 

3.  Compare  the  performance  of  the  constant-area,  supersonic- 
supersonic  ejector  with  that  of  the  constant-area,  subsonic- 
super  '.onic  ejector  as  applied  to  high-energy,  chemical  laser 
systems . 


Diluent 


(H  F ~ 2. 7 it,  DF  ~ 3.8  u)  1 ^owftr  Atf?)os< 

Figure.  1.0-1  Typical  Flow  n;agra,„  for  a m,,,.,,  _ 
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Figure  1.0-2  Laser  Cavity  Scf 
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Figure  1.0-3  Subsonic-Supersonic  Pumping  System 


Figure  1.0-4  Supersonic-Supersonic  Pumping  Sys'er.: 
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2.0  THEORETICAL  ANALYSIS  OF  THE  CONSTANT- AREA, 

SUPERSON IC - SUPERSON I C EJ ECTOR 

In  the  constant-area,  supersonic -supersonic  ejector  of  Fig.  2.0-1, 
the  primary  and  secondary  streams  enter  at  supersonic  Mach  numbers  where 
they  begin  to  interact,  mix,  and  diffuse.  Proceeding  through  the  ejector, 
thc%e  streams  continue  to  mix  and  diffuse,  thus  approaching  a uniform 
flow  at  the  ejector  exit.  The  resuit  of  this  interaction,  mixing,  and 
diffusion  process  can  be  described  by  a one-dimensional,  compressible 
flow  model  consisting  of  two  components,  an  overall  analysis  of  the 
constant-area  mixing  section,  stations  1 to  3,  and  an  analysis  of  the 
nearly  inviscid  interaction  region  just  downstream  of  the  confluence  of 
the  primary  and  secondary  streams. 

By  tradition,  the  governing  equations  fo~  subsonic-supersonic  ejector 
models  are  nondimensional ized  by  the  appropriate  primary  flow  variables 
to  allow  for  a zero  secondary  mass  flow  rate,  and  this  tradition  is  main- 
tained in  the  following  analysis,  however,  in  preference  to  people  in  the 
chemical  laser  field,  the  results  arc  presented  as  nondimens iona 1 i zed  by 
the  secondary  flow  variables. 

2.1  ONE-DIMENSIONAL  OVERALL  MIXING  SECTION  ANALYSIS 

By  hypothesis,  the  streams  entering  the  ejector  must  be  supersonic. 
This  criteria  restricts  the  ejector  operation  to  a plane  as  illustrated 
»n  Fig-  2.1-1,  and  prescribes  the  boundaries  of  this  plane.  So  long  as 
the  entering  streams  remain  supersonic,  the  ejector  mass-flow  ratio  W /W 
is  established  independent  of  the  ejector  exit-plane  pressure  but  directly 
proportional  to  ti.e  primary-to-sccondary  static  pressure  ratio  Ppi/Rsi  at 
the  confluence  point.  The  right -most  ooundary  of  the  plane  of  operation 
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defines  the  maximum  compression  ratio  and  is  determined  by  applying  the 
:onservation  equations  to  the  control  volume  of  Fig.  2.1-2  together  with 
the  following  assumptions: 

(1)  Steady  flow,  2-  H 0. 

ot 

(2)  Piecewise  uniform  flows  at  station  1 and  uniform  flow  at 
station  3. 

(3)  The  primary  and  secondary  gases  obey  the  perfect  gas 
relationships. 

(4)  The  primary  and  secondary  streams  mix  ideally  to  form  a mixed 
gas  at  station  3. 

(5)  Negligible  wall  shear  stresses. 

(6)  Adiabatic  flow  between  stations  1 and  3. 

(7)  Nc  shaft  or  shear  work  between  sl^tions  1 and  3. 

(8)  Negligible  body  forces. 

(9)  The  flow  in  the  primary  nozzle  is  isentropic  from  its 
stagnation  state  to  the  state  at  station  1. 

The  fundamental  equations  of  continuity,  momentum,  in  the  flow 
direction,  and  energy  ?re,  respectively: 


cv  cs 


I F*  “ 


3 

3t 


^Jvx(  pdv)  +j>' 
cv  cs 


Vx(Pv 


dA) 


(2.1-2) 
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0Q  0Wss  a r V2 

W ' W~  = IT J (u  + T + gz 


’ ) (pdV) 


+ f(h  + ~ * 


gz  + • ••) (pV  • dA) 


(2.1-3) 


Applying  assumptions  (1,2)  to  Eqn.  (2.1-1)  yields 


ps  i Vs  i As  l ' Pp  i Vp  l i + ^f.c  VM3  \a  0 


or  in  terms  of  the  mass  flow  rate  W = pAV, 


WM  WC 

M , S 

\u  +wT  • 


Applying  assumptions  (1,2,5)  to  Eqn.  (2.1-2)  gives 


PA  + P A_  - P A = - V (0  V A ) 
S!  si  PI  >1  M3^13  si  lMsi  si  SlJ 


- V (p  V A ) + V (p  V A ) , 
PI  VMP1  P1P1J  M3  [.C  IB'  * 


or  in  terms  of  the  Mach  number, 


ps,Asi»  * Vs,>  * PpA.o  - V$.>  ■ * v4> 


Since  A^  + A^  = A^  for  a constant-area  mixing  tube,  the  result  is 


P P 

M3  P 1 


Apj  * *-^S  ’^l  ^ + fi*-Yp,f>Vi^ 


1 + A^-J  MW 


where  fj  (y,M)  = 1 + yM2  . 


(2.1-6) 
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Equation  (2.1-3)  together  with  assumptions  (1,2,6, 7,8)  yields 


^hsi  + 2 ^PsiVsiAsi^  " (hri  + 2 -^ppi  Vpi 


+ 2 ■*  VM3  Am3  " ° 


or  in  terms  of  the  mass  flow  rate  and  stagnation  enthalpy  hQ  = h + — , 


h W + h K - h W 
sos  pop 


But  since  ho  = CpTo  + constant  for  an  ideal  gas,  the  result  is 


Tmo  Wm  «^m  Tso  Ws  (CP)S 

Tpo  * Wp  ' (Cp)p  ' + Tpo  * Wp  ' (Cp)p 


Using  assumptions  (3,4),  the  specific  heats  of  the  mixed  flow  are 
related  to  their  primary  and  secondary  stream  counterparts  by 


Vcp)M  = Ws  + Wr  » 


Wm  = Ws  + Wp  • 


Rearranging  Eqn.  (2,1-8)  as 


WM  (Cp)«  WS 

WP  ‘ «^P  _ + WP  ' ’ 


v p 

and  noting  that  Cp  = (^y)  ^or  a rer^ect  gas,  gives  the  useful 
relations 
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(Vs  _ Ys  rYp  - 

«Vp  ~ M«s  1ys  'Ml  Yp  J 


. _<Vm  \ < ^P_  f YS  fYP  - 1 

Wp  (Vr  ’ + WP  ’ Mws  IYS  - M l Yp 


Similarly,  rearranging  Eqn.  (2.1-9)  as 


"m  «Vm  , *s  <Cv>s 
\ ' “p  ■ (cvY 


1 ^ R 


and  noting  that  Cy  = (— y)  ^ for  a perfect  gas,  gives  the  useful  rela- 


tions 


(Vs  Mwp  ft  * 1 


(CV)P  Mwg  Iys  - 1J  ’ 


(2.1-12) 


WM  (Cv>m  , Ws  MwP  fYp  * M 

Wp  (Vp  " + Wp  ’ Mws  Iys  - l) 


P s '•'s 


(2.1-13) 


The  mass  flow  rate,  W,  is  expressed  in  terms  of  the  mass  flow  func- 


tion by 


TK  ^ * To]V2  = M {y[1  + y^2  = f2(Y»M)  • (2.1-14) 


Then  using  this  relation,  the  secondary-to-primary  mass  flow  ratio  is 
related  to  the  static  pressure  ratio  at  the  confluence  point  by 


Ws  PS1  ^ \l  rMWS  ^ V..V2  ^2  (YS  1 ^ 

% V \l  YS0  M (yp  ’^P 1 ^ 
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Since  y = for  a perfect  gas,  the  mixed  gas  property  is  given  by 


(CP)M  (Cp)p 

Ym  " " (CA 


(Cr)M 

(Cp)F 

CVr 


which,  with  Eqns . (2.1-11)  and  (2.1-13),  becomes 


'm  W. 


W„  Mw„  y„ 

_s_  . , s N 'p  ' 

W Mw_  To  - 1 j + - 1 

PS'S  'P 

Mw  Y Y 

p * C 1 p 

*7  [(vr>  - 1]  * [(v^> 


Subtracting  Eqn.  (2.1-8)  from  Eqn.  (2.1-9)  and  using  the  perfect  gas 
relation  — = Cp  - Cy  yields 


IV  Mw  W Mw 

M P S P 

>V  Mw  + W Mw 

P M V S 


Then  applying  Eqn.  (2.1-4)  and  rearranging  results  in 


^+1 


Mw 

W. 

Mw_ 

' 

p 

b 

, + i 

Mw 

W 

p 

s 

Substituting 

relations  (2.1 

-10) 

and 

(2- 

T 

W 

Mw_ 

Y 

1 

so 

r 

p 

’s 

] 

‘r»J0 

T 

PO 

" t * 

Mws 

U • 

lj 

T~~ 

W„  Mw 

( 

Y v 

po 

s p 

. - r.  • ■ 

S 

Wp  Mws 

-Ys 

- lJ 

T 

IYp  - 1] 


(2.1-18) 
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Using  the  mass  flow  function  (2.1-14),  the  mixed-to-primary  mass  flow 
ratio  is  expressed  as 


_ PM3  \a  MWfC  ^PO  j/2  ^2  (YM’^M3  ^ 

“p  ' pp.  ' \7  tiS7 ' ^3  pjv  • 


Applying  Eqn.  (2.1-4)  and  the  relation  + Ap  = A^n  , this  equation 
becomes 


Pm3  ri  + ^5  1 ■,  r^M  ^P0nl/2 

Pn  S’V  ^V^TT 


, WS 
+ \ 


Then  eliminating  p^.  with  Eqn.  (2.1-5)  yields 
hr  l 


fjW- 


Mw  T PA 

r M *P0  1/2  .SI  'si 

Cf 5T  * ~ 1 rp~  * t~  * fs  CYa-^i)  + fi  (V  NVi)] 

P 7.B *P  1 >1  1 r 1 

% 

t 1 + f2  *-Yp  ,mpi  ) 

P 


(2.1-19) 


where 


r3  (Y,M) 


f. (V-M) . i.y 


f r 


2 (Y > M)  M \ Y r 1 + (^I)M2  ]} 


(2.1-20) 


The  exit  Mach  number  is  obtained  by  solving  the  quadratic  equation 


(Y„-D 


C~T ^ - ym]  (ML)2  + [f^  - ~ ^ = 0 • (2-1-21) 


Tnen  the  mixed-to-primary  static  pressure  ratio  is  found  by  solving, 
in  order: 


Wp  ’ 


(1)  Equation  (2.1-15)  for 
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(2)  liquation  (2.1-16)  for  ym  » 


Mw. 


(3)  Equation  (2.1-17)  for 


M 


Mwp  ’ 


(4)  Equation  (2.1-18)  for  ^ 


MO 


PO 


(5)  Equation  (2.1-19)  for 


(6)  Equation  (2.1-21)  for  Mf>c  , and 


(7)  Equation  (2.1-5)  for 


M3 


PI 


Mws  Tso 


S 1 


where  y$  , yp  , ~ — , Msl  , , j-— , and  ^ — are  the  independent 

1 p PO  ‘ PI  PI 

variables  or  ejector  parameters.  For  computational  purposes,  the  static 
pressure  Pp,  may  be  obtained  from  the  stagnation  pressure  Pp0  using 
assumption  (9)  and  the  isentropic  flow  relation 


~ (y,M)  = [l  + (^-)  M2]7'1 


(2.1-22) 


For  chemical  laser  applications,  relation  (2.1-22)  may  not  be  applied  to 
the  secondary  stream  unless  P$0  is  taken  at  station  1 in  Fig.  2.0-1  since 
the  laser  cavity  flow  will  not  be  isentropic. 


2 

It  should  be  noted  that  Eqn.  (2.1-21)  has  two  roots  for  giving 
subsonic  and  supersonic  values  of  . The  right-most  boundary  of  the 
plane  of  operation.  Fig.  2.1-1,  is  calculated  from  the  subsonic  value  of 
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diffused  through  a normal  shock  wave,  they  yield  exactly  the  subsonic 
solutions  along  the  right-most  boundary  of  the  plane.  Of  course,  the 
ejector  may  operate  anywhere  within  the  boundaries  of  the  plane  of 
supersonic-supersonic  operation;  howev  r,  given  Wp/Ws  or  PF0/psl  > the  most 
desirable  operation  is  at  the  right-most  boundary  since  the  potential 
exists  for  operation  at  this,  the  maximum  compression  ratio. 

2.2  ONE-DIMENSIONAL  ANALYSIS  OF  THE  INVISCID  INTERACTION  REGION 

The  upper  boundary  of  the  plane  of  operation,  Fig.  2.1-1,  is  also 
dictated  by  the  requirement  that  both  the  secondary  and  primary  streams 
remain  supersonic.  If  the  primary-to-secondary  static  pressure  ratio  at 
the  confluence  point  is  greater  than  1 (P  /P  >1),  the  secondary  flow 
is  compressed  by  the  mutual  interaction  of  the  primary  and  secondary 
streams  within  the  mixing  tube.  This  process  is  limited,  in  a one- 
dimensional sense,  to  a "nearly"  reversible  recompression  to  sonic  flow 
at  the  minimum  area  as  determined  by  the  control  volumes  shown  in  Fig. 

2.2-1.  Thus,  the  constant-area,  supersonic-supersonic  ejector  couples 
the  effect  of  an  ideal  aerodynamic,  supersonic  diffuser  and  momentum 
transfer  through  viscous  mixing. 

The  control  volume  of  Fig.  2.2-1 (a)  extends  from  station  1 to 
station  2.  In  addition  to  the  assumptions  listed  in  Section  2.1,  the 
following  additional  assumptions  are  made: 

(10)  The  streams  remain  distinct  and  do  not  mix  between  stations 
1 and  2 . 

(11)  The  flow  is  isentropic  for  each  stream  between  stations  1 and  2. 

(12)  The  average  pressures  of  the  streams  can  be  different  at  each 
cross-section;  thus,  continuity  of  static  pressure  at  the 
boundary  between  the  streams  is  not  satisfied  by  this  flow  model. 
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(lo)  The  Me eh  number  of  the  secondary  flow  at  station  2 is  M^2  - 1 • 
(14)  The  static  pressures  are  such  that  P > P$i . 

For  an  isentropic,  compressible  flow,  the  area  ratio,  A/A*,  is 
expressed  in  terms  of  the  area  ratio  function  by 


A_ 

A* 


* <^>  f4(Y,„, 


(2.2-1) 


since  Mg2  = 1,  Ag2  = A*  and  for  a constant-area  mixing  tube 

AS1  + \l  = AS2  + \2-  ^ 


\2  Asi  As* 

{1  + _1L  [1  . _£_]}  , or 


\*  \*  Api 


S 1 


Ap 


A*  =f4(YP,NVi}  {1  + Apj  n ' (f4(Ys,Msl))])  ; (2‘2*2) 


and 


Ap2 

f4  *■  Yp  ’ 2 ^ ~ A*- 


(2.2-3) 


can  be  solved  for  the  supersonic  value  of  . 

By  assumptions  (10,11),  W , Wp , T$0  , and  Tp0  are  constant  from 
station  1 to  station  2 in  Fig.  2.2-1 (a).  Then  the  mass  flow  function 
(2.1-14)  gives 


S 2 


SI 


S 2 


SI 


f2<VMS,> 

f2(VMS2> 


and 


(2.2-4) 


Pp  2 # ^2  _ f2  *■  YP  ’ 1 ^ 

Ppi”  AP1  = f2(V^P 


(2.2-5) 
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Applying  the  momentum  equation  (2.1-2)  to  the  combined  control  volume 
of  Fig.  2.2-l(b)  together  with  assumptions  (1,2,5,10)  yields 

PA  + P A -PA  _pa_  = - V (pVA„) 

SI  SI  PI  PI  S 2 S 2 P2P2  SI  lMSl  SI  Si' 


- V (p  V A ) + V (p  V A ) + V (p  V A ,) 

PI  VKP1  P1P1-'  S2  ^S2  S2  S2-'  P2  VIP2  P2  P2  ' 


or  in  terms  of  the  Mach  number. 


ps,\,»  * v£.>  * ■VAP1  * = 


psA2<>  * V&>  * pP2\2(i  ♦ yX2) 


Then  using  the  function  (2.1-6),  the  result  is 


PS1  ^ AS1  fl  ^YS  ,MS1  + PS2  _ AS2  > PS1  AS1  fl  ^YS  ,MS2  -* 

Pp  l \ l f l ’ % i ^ Ps  l As  i Pp  1 \ l fi  l •* 


PP2  2 ^1^P’^P2^ 

^7  ■ \T f! 


(2.2-6) 


Combining  Eons.  (2.2-4),  (2.2-5),  and  (2.2-6)  with  = 1 yields 
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where  f (Y,M)  is  the  function  (2.1-20). 

Then  the  static  pressure  ratio  Psl/Ppi  for  an  isentropic  recompres- 
sion of  the  secondary  stream  to  sonic  conditions  at  station  2 is  obtained 
by  solving,  in  order: 
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(1)  Equation 

(2)  Equation 

(3)  Equation 

where  yg , yp , MS1 , 

parameters;  and  P 
function  (2.1-22). 
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(2.2-2)  for  , 

(2.2-1)  for  Mp2  > 1 , and 


(2.2-7)  for 


Mpj  , and  — — arc  the  independent  variables  or  ejector 
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may  again  be  related  to  Ppo  by  the  isentropic  fl 


ow 


2.3  PARAMETRIC  RESULTS 

The  maximum  compression  ratio  for  a given  supersonic-supersonic 
ejector  configuration,  i.e.  given  Yg , yp  , Mw$/Mwp,  Tgo/Tpo,  Mgl  , Mpj  , and 
A‘j l / Ap  i > *s  defined  by  the  intersection  of  the  upper  boundary  or  "upper 
limit  line"  of  the  plane  of  operation  as  calculated  from  Section  2.2,  and 
the  right-most  boundary  of  the  plane  as  calculated  from  Section  2.1.  This 
intersection  point  is  termed  the  "upper  limit  point"  or  "ULP"  in  Fig. 

2.1-1  . 


The  static  pressure  ratio  P /P  may  take  on  values  less  than  1; 
however,  the  roles  of  driver  and  driven  streams  are  interchamged  and  the 
plane  of  operation.  Fig.  2.1-1,  would  be  reproduced  with  the  primary  and 
secondary  subscripts  interchanged.  Then  for  practical  purposes,  the  lower 
boundary  of  the  plane  of  supersonic-supersonic  operation  was  taken  to  be 
the  "marched  pressure  line"  where  Ppi  = P . The  intersection  of  the 
matched  pressure  line  with  the  right-most  boundary  of  the  plane  is  termed 
the  "matched  pressure  point"  or  "MPP"  in  Fig.  2.1-1. 

The  influence  of  the  input  variables,  yg  , yp  , Mwg/Mwp,  Tgo/Tp(), 

M , Nip j , and  ^S1/\l  > on  the  plane  of  operation  is  illustrated  in 
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Figs.  2. 3-1  to  2.3-6.  The  procedure  for  producing  each  figure  was  to  hold 
all  the  ejector  parameters  at  constant  values,  save  one,  and  to  vary  this 
parameter  over  a wide  range.  Rather  than  plot  the  resultant  planes  of 
operation,  only  the  loci  of  upper  limit  and  matched  pressure  points  are 
given.  The  base  values  for  the  input  variables  were  in  all  figures: 

Ys  - i.40,  Mws/MWp  = 1.00,  yp  = 1.40,  T so/Tpo  = 1.00, 

Msi  = 2‘00’  AS1/\1  = 2-00’  Mp!  = 4-00> 

while  each  point  on  the  graphs  represent  valid  solutions  to  the  equations 

of  Sections  2.1  and  2.2,  the  figures  with  y , y , and  Mw  /Mw  as  the 
abcissa  may  not  represent  physically  realistic  solutions  except  at  the 
base  value,  indicated  by  the  vertical  dashed  line,  since  only  limited 
combinations  of  y and  Mw  occur  in  nature.  It  should  also  be  noted  that 
the  identical  influence  of  Mv/S/Mwp  and  Tso/Tpo  on  the  plane  of  operation, 
as  shown  in  Fig.  2.3-3,  is  simply  a coincidence  of  the  choice  of  base 
values.  The  most  significant  information  given  by  Figs.  2.3-1  to  2.3-6 
is  the  slope  of  each  curve  at  the  base  value  or  vertical  dashed  line,  since 
this  slope  indicates  the  partial  derivative  of  the  dependent  variable, 
Ppo/Psl,  Wp/Wg,  or  P:,rr/PSi  > with  respect  to  the  independent  variable. 

2.4  COMPUTER  PROGRAMS 

The  analyses  of  Sections  2.1  and  2.2  were  the  bases  for  the  develop- 
ment of  computer  programs  for  analyzing  the  performance  of  constant-area, 
supersonic-supersonic  ejectors.  These  programs  along  with  sample  input 
and  output  are  presented  in  detail  in  APPENDICES  7.2  and  7.3  for  programs 
CASSE  and  CASSEP,  respectively. 

These  computer  programs  were  used  to  calculate  the  parametric  results 


presented  in  Section  2.3. 
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Figure  2.3-3  Influence  of  Mw  /Mw 
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3.0  EXPERIMENTAL  ANALYSIS  OF  THE  CONSTANT-AREA, 
SUPERSONIC-SUPERSONIC  EJECTOR 

A series  of  small-scale,  cold-flow,  air-to-air  ejector  experiments 
were  conducted  to  verify  the  theory  of  Section  2.0  and  to  provide  a more 
detailed  description  of  the  flow  field  than  necessarily  allowed  by  a one- 
dimensional analysis. 

3.1  EXPERIMENTAL  APPARATUS 

A half-section  of  the  axisymmetric , supersonic-supersonic,  ejector 
model  is  given  in  Fig.  3.1-1.  This  ejector  model  was  designed  and  fabri- 
cated to  facilitate  a rapid  replacement  of  nozzles  and  mixing  tubes. 

Primary  flow  enters  the  ejector  from  a large  stagnation  chamber  and  is 
accelerated  through  an  elliptical  entrance  section  into  the  interchange- 
able nozzle.  Secondary  flow  enters  the  secondary  stagnation  chamber  from 
two  sides  and  is  accelerated  between  the  nozzle  wall  and  elliptical 
entrance  section  at  the  mixing  tube  base.  Static  pressure  taps  were  in- 
stalled along  the  mixing  tube  wall  from  the  nozzle  exit  plane  to  the 
mixing  tube  exit  at  increments  of  one  tube  radius.  A series  of  six  static 
pressure  taps,  not  shown,  were  equally  spaced  about  the  mixing  tube  axis 
at  the  nozzle  exit  plane  to  check  the  concentricity  of  the  nozzle  and 
mixing  tube.  In  addition,  five  static  pressure  wall  taps,  not  shown,  were 
located  upstream  of  the  nozzle  exit  plane  at  increments  of  2.54  mm  to 
ensure  that  the  entering  secondary  flow  was  supersonic. 

For  experimental  purposes,  the  supersonic  primary  and  secondary  flows 
were  produced  by  concentric  nozzles  of  the  continuous-slope  type,  as  shown 
m Fig.  3.1-2.  The  primary  stream  was  directed  through  the  central  nozzle, 
while  the  secondary  flow  was  accelerated  through  the  annular  passage 
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formed  by  the  outer  nozzle  wall  and  the  constant-area  mixing  tube.  Nozzle 
and  mixing  tube  specifications  are  also  given  in  Fig . 3.1-2  with  the  cor- 
responding dimensionless  parameters  summarized  in  Table  3.1-1.  While 
confined  to  an  axisymmetric  geometry  and  narrow  Mach  number  range  due  to 
test  facility  limitations,  these  experimental  configurations  should  be 
adequate  for  an  initial  evaluation  of  supersonic-supersonic  ejector  per- 
formance . 

The  continuous-slope  nozzles  constructed  for  these  experiments  were 
designed  by  the  method-of-characteristics  to  produce  uniform  velocity  and 
pressure  distributions  at  the  confluence  of  the  secondary  and  primary 
streams,  thus  satisfying,  at  least  ideally,  the  assumptions  for  the  one- 
dimensional theory  of  Section  2.0.  A photographic  enlargement  of  one  of 
these  nozzles  emphasizing  the  secondary  or  outer  nozzle  wall  profile  is 
presented  in  Fig.  3.1-3.  As  a design  simplification,  all  four  nozzles 
were  machined  to  the  same  length  which  necessitated  the  addition  of 
constant-area  sections  at  the  nozzle  throat  and  exit  segments  as  specified 

Table  3.1-1  Dimensionless  Parameters  for  the  Experimental, 
Supersonic-Supersonic  Ejector  Model 
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in  Fig.  3.1-4.  Although  these  constant- area  sections  would  have  no  effect 
for  a truly  inviscid  flow,  they  would  promote  boundary  layer  growth  under 
actual  experimental  conditions.  Figure  3.1-5  is  a near  actual  size 
photograph  of  all  four  supersonic-supersonic  nozzles  showing  their  simi- 
larity in  size  and  shape. 

Photographs  of  the  ejector  components  are  given  in  Fig.  3.1-6  showing 
both  front  and  rear  views  of  the  secondary  stagnation  chamber.  It  should 
be  noted  from  Fig.  3,1-2  that  the  two  constant-area  mixing  tubes  are  10 
tube  diameters  in  length  and  that  the  interchangeable  nozzles  were  con- 
structed with  a separate  elliptical  entrance  section.  A view  of  the 
ejector  model  in  a partially  assembled  state  is  given  in  Fig.  3.1-7 
showing  the  position  of  the  supersonic-supersonic  nozzle  in  the  secondary 
stagnation  chamber,  one  of  two  side  entrance  ports  for  the  secondary  flow, 
and  one  of  two  "U”-shaped  baffles,  the  visible  one  surrounding  the  secon- 
dary stagnation  pressure  tap,  which  prevent  the  formation  of  a large 
circular  vortex  about  the  nozzle  axis.  All  the  ejector  components  were 
precision  machined  to  design  specifications  with  less  than  1.27  mm 
clearance  and  0-ring  seals  at  all  the  mating  surfaces  to  ensure  proper 
alignment  and  sealing. 

Figure  3.1-8  is  a flow  diagram  for  the  experimental  program.  Air 
from  a common  supply  branched  to  two  automatic  control  valves  which  main- 
tained a constant  stagnation  pressure  in  the  primary  and  secondary 
chambers.  The  air  then  passed  through  standard  VDI  nozzles  for  measure- 
ment of  the  primary  and  secondary  mass  flow  rates.  A final  stagnation 
chamber  was  located  at  the  mixing  tube  exit  immediately  upstream  of  a 
back  pressure  control  valve  to  prevent  local  disturbances  at  this  valve 
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from  influencing  the  velocity  and  pressure  distributions  of  the  mixed, 
subsonic  flow  at  the  tube  exit.  The  back  pressure  control  valve  was  a 
sliding  block  arrangement  with  two  blocks  closing  symmetrically  across 
the  exit  duct. 

Figure  3.1-9  is  a photograph  of  the  small-scale,  ejector  model 
installed  on  the  test  chamber,  which  served  as  the  primary  stagnation 
chamber,  with  a balance  handle  for  manual  adjustment  of  the  back  pressure 
control  valve  visible  in  the  upper  portion  of  the  photo.  Also  visible 
are  the  pressure  lines  leading  to  the  mixing  tube  wall  taps  and  a portion 
of  the  silencer  which  supports  the  final  stagnation  chamber  and  back 
pressure  control  valve.  For  one  experiment,  a traversing  pitot  probe, 
as  shown  in  Fig.  3.1-10,  was  added  between  the  mixing  tube  and  final 
stagnation  chamber  for  measurement  of  the  exit  Mach  number  distribution. 

3.2  EXPERIMENTAL  PROCEDURE 

The  experimental,  supersonic-supersonic  ejector  model  was  installed 
on  a test  chamber  of  the  continuous  flow  facility  in  the  Mechanical 
Engineering  Laboratory  and  was  operated  with  dry,  compressed  air.  The 
tests  were  run  at  secondary  stagnation  pressures  of  98  to  270  kPa  abso- 
lute, primary  stagnation  pressures  of  269  to  741  kPa  absolute,  and 
stagnation  temperatures  of  approximately  294  K.  The  primary  nozzle 

5 

Reynolds  number,  based  on  the  throat  diameter,  varied  from  4.3  x 10 
at  PpQ  = 269  kPa  to  1.2  x 10&  at  Fpo  = 741  kPa.  Secondary  nozzle  throat 
Reynolds  numbers  are  indicated  in  the  experimental  results. 

Pressure  data  were  taken  from  Bourdon-tube  gauges  and  manometers 
or  with  a strain  gauge  transducer-digital  counter  system.  The  primary 
and  secondary  stagnation  pressures  and  static  pressures  upstream  of  the 
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standard  VDI  nozzles  were  measured  with  Bourdon-tube  gauges  except  at 
levels  near  atmospheric  pressure,  in  which  case  a mercury  manometer  was 
used.  Pressure  differences  across  the  standard  VDI  nozzles  were  measured 
with  either  mercury,  Meriam  3*",  or  water  filled  U-tube  manometers,  depen- 
ding on  the  magnitude  of  the  pressure  difference.  Photographic  records 
of  all  static  pressures  along  the  mixing  tube  wall  were  taken  from  mercury 
manometer  boards.  Stagnation  pressure  readings  lor  the  traversing  pitot 
probe  were  taken  with  the  strain  gauge  transducer-digital  counter  system. 

Although  six  static  pressure  taps  were  available  for  the  measurement 
of  P , only  one  was  used  in  each  experiment,  and  this  pressure  was  read 
with  both  a mercury  manometer  and  a strain  gauge  transducer.  An  attempt 
was  made  to  measure  an  average  value  of  P by  manifolding  all  six  pres- 
sure  taps;  however,  the  resultant  readings  were  lower  than  for  any  one 
given  pressure  tap. 

Since  the  stagnation  temperatures  Tpo  and  Tg0  are  primarily  dependent 
on  the  supply  temperature,  the  primary  stagnation  temperature  was  measured 
with  a dial  thermometer  and  values  for  Tg0  and  the  temperature  of  the  air 
upstream  of  the  standard  VDI  nozzles  were  assumed  to  equal  this  value  for 
T . 

PO 

Maximum  compression  ratio  data  for  each  experimental  ejector  con- 
figuration were  obtained  in  the  following  manner.  With  the  back  pressure 
control  valve  fully  open,  the  secondary  stagnation  pressure  P$0  was  set 
at  a fixed  value  to  be  maintained  by  the  secondary  automatic  controller 
(see  Fig.  3.1-8).  The  automatic  controller  for  the  primary  stagnation 

^Product  of  the  Meriam  Instrument  Company,  Cleveland,  Ohio. 
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pressure  PpQ  was  then  , t at  a fixed  value  such  that  the  static  pressure 

ratio  P /P  would  lie,  at  least  theoretically,  in  the  range  between  the 

matched  pressure  and  upper  limit  lines.  The  back  pressure  valve  was  then 

closed  until  the  value  of  P , as  indicated  by  the  digital  counter,  began 

to  rise,  at  which  point  the  mercury  manometer  board  was  photographed. 

This  process  was  repeated  for  different  combinations  of  Pg0  and  PpQ  until 

the  full  range  of  P /P  was  covered  unless  otherwise  restrict^'.  by  the 
pi  si 

maximum  supply  pressure. 

While  the  above  process  is  necessarily  transient  in  nature,  the 
results  are  thought  to  be  good  since  the  point  at  which  P begins  to 
rise  rapidly,  indicating  a transition  from  supersonic  to  subsonic  flow, 
was  quite  well  defined  and  the  time  constant  for  the  mercury  manometer 
boards  is  much  greater  than  for  the  strain  gauge  transducer-digital 
counter  system.  Thus,  measurements  of  P_,  and  P,_  from  the  mercury  man- 
ometers  should  be  very  close  to  the  limit  of  supersonic-supersonic 
operation.  This  thought  is  a] sc  supported  on  the  basis  of  data  repeat- 
ability which  was  about  +_  1%  for  Prc/Psi  with  = 1.50. 

Traverses  with  the  pitot  probe  at  the  mixing  tube  exit  could  not  be 
completed  at  exactly  maximum  back  pressure  conditions.  The  above  proce- 
dure was  followed  until  P began  to  rise,  and  then  the  back  pressure 
valve  was  opened  slightly  to  reach  a stable  state  for  the  traverse. 

3.3  EXPERIMENTAL  RESULTS 

Maximum  compression  characteristics  for  the  experimental,  constant- 
area,  supersonic-supersonic  ejector  configurations  listed  in  Table  3.1-1 
are  presented  in  Figs.  3.3-1  to  3.3-4.  In  each  figure,  the  data  points 
are  plotted  on  the  plane  of  supersonic-supersonic  operation  as  predicted 
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by  the  one-dimensional  analysis  of  Section  2.0  and  illustrated  m Fig. 
2.1-1.  These  data  points  indicate  the  maximum  compression  ratio,  P^/Pgj  > 
for  subsonic  values  of  and,  therefore,  should  lie  on  the  right-most 
boundary  of  the  plane  of  supersonic-supersonic  operation.  Also,  given 
the  compression  ratio  Prc/Psl , the  mass  flow  ratio  data  points  should  lie 
on  the  pressure  ratio  data  points,  i.e.  the  triangles  should  lie  on  the 
circles  directly  below  them,  since  Wp/Ws  is  theoretically  proportional  to 
Ppo/Psl.  In  all  cases,  given  Ppo/Psl»  the  experimental  values  for  the 
maximum  compression  ratio  P^j/Pgj  were  15  to  18  percent  less  than  the 
theoretical,  except  for  = 1.50,  in  which  case  the  error  was  21  to  22 
percent,  and  the  experimental  values  of  Wp/W$  were  0 to  44  percent  greater 
than  the  theoretical.  It  should  also  be  noted  that  the  percent  errors  in 
PM3/P$1  and  especially  Wp/Wg  increase  with  a decrease  in  the  secondary 
nozzle  throat  Reynolds  number  as  defined  by 
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Wall  pressure  distributions  at  maximum  compression  conditions  for 
the  ejector  configurations  listed  in  Table  3.1-1  are  given  in  Figs.  3.3-5 
to  3.3-8.  Three  pressure  distributions  are  plotted  in  each  figure  which 
correspond  to  data  points  near  the  matched  pressure  point,  upper  limit 
point,  and  at  some  point  midway  between  in  Figs.  3.3-1  to  3.3-4. 

At  the  matched  pressure  point,  the  wall  pressure  is  initially  con- 
stant. followed  by  a near  linear  rise  whi<_h  levels  off  at  the  mixing  tube 
exit.  The  linear  rise  covers  a smaller  portion  of  the  nixing  tube  at  the 
larger  values  of  M corresponding  to  the  higher  initial  velocity  of  the 
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secondary  stream  and  the  fact  that  viscous  mixing  of  the  primary  and 
secondary  streams  is  accomplished  in  a shorter  length  as  Mg ^ approaches 

NVi- 

At  the  upper  limit  point,  the  secondary  stream  undergoes  a large 
initial  recompression  followed  by  a near  linear  rise  which  again  levels 
off  at  the  mixing  tube  exit.  The  initial  pressure  rise  can  be  attributed 
tc  the  aerodynamic  diffuser  formed  by  the  primary  stream  as  it  expands 
into  the  mixing  tube.  As  noted  for  the  matched  pressure  data,  the  linear 
portion  of  the  pressure  rise  occupies  a smaller  segment  of  the  mixing  tube 
as  ML  approaches  Mpj  , though  the  total  pressure  rise  is  much  greater 
since  more  energy  is  transferred  to  the  secondary  stream  as  Ppo/psl  is 
increased . 

The  wall  pressure  distributions  of  Figs.  3.1-5  to  3.1-8  also  show 
that  the  low  experimental  values  for  the  maximum  compression  ratio  in 
Figs,  3.1-1  to  3.1-4  may  be  due,  in  part,  to  the  mixing  tube  length 
which  was  equal  to  10  tube  diameters  in  nil  the  experimental  ejector  con- 
figurations, 8 to  12  tube  diameters  being  sufficient  for  most  subsonic- 
supersonic  ejectors.  If  the  mixing  cube  is  long  enough  for  completion  of 
the  viscous  mixing  process,  yet  not  too  long  such  that  wall  friction 
becomes  important,  then  the  wall  pressure  distributions  should  exhibit  a 
near  linear  rise  which  levels  off  to  i constant  value  as  the  primary  and 
secondary  streams  mix  to  a uniform  flow.  With  the  exception  of  Fig.  3.3-8, 
where  ^ = Ni,  , it  appears  that  the  linear  portion  of  the  wall  pressure 
distributions  end  at  the  mixing  tube  exit,  that  the  viscous  mixing  and 
diffusion  process  may  not  have  been  completed,  and  that  the  compression 
ratio  might  oe  increased  with  a longer  mixing  tube.  This  conclusion  is, 
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moreover,  supported  by  the  exit  Mach  number  distributions  at  maximum  com- 
pression conditions  in  Fig.  3.3-9  which  show  that  for  M = 2.00,  the 
potential  exists  for  further  mixing  and  diffusion  even  for  Ppi/Pgl  near 
the  matched  pressure  line. 

The  large  percent  errors  in  Wp/W  and  the  changes  in  compression  char- 
acteristics with  ke  noted  in  Figs.  3.3-1  to  3.3-4  are  probably  due  to  sep- 
aration of  the  secondary  flow  upstream  of  the  confluence  point.  In  the  one- 
dimensional  theory  of  Section  2.0,  it  is  assumed  that  the  secondary  stream  would 
undergo  an  isentropic  recompression  to,  at  the  least,  sonic  conditions  for  P <_ 

P as  shown  schematically  in  Fig.  3.3-10(a).  In  the  real  situation,  however,  the 
secondary  stream  may  not  be  able  to  change  directions  at  station  1 without  separa- 
ting the  boundary  layer  as  illustrated  in  Fig.  3.3-10(b).  The  flow  losses  through 
the  oblique  shock  structure  may  not  be  significant  even  with  secondary  separation, 
but  the  irea  is  effectively  reduced  with  a corresponding  increase  in  Pgi  and 
decrease  in  M$1  ; hence,  larger  secondary  stagnation  pressures  and  mass  flow  rates 
are  required  to  obtain  the  desired  values  of  Pp  /Pgi  . The  region  of  secondary  sep- 
aration should  grow  with  increasing  values  of  Ppj  /Pgi  and.  decreasing  secondary  Rey- 
nolds number,  ana  this  fact  is  evident  in  Fig.  3.3-11  where  M . as  calculated  from 

° si 

p a 

the  isentropic  flow  function  p (y,M),  is  plotted  versus  / Pg ; for  the  experimer 
tal  data  of  Figs.  3.3-1  to  3.3-4.  Some  methods,  such  as  Zukoski's  ( lb) are  avail- 
able for  description  of  the  separation  region,  but  prediction  of  the  static  pres- 
sure ratio  P /P  at  the  actual  onset  of  separation  requires  empirical  data, 
thus  defeating  most  simple  models  for  parametric  evaluation  of  supersonic-super- 
sonic ejector  performance. 

•j* 

The  results  of  an  experimental  investigation  of  a typical  ejector  config- 
uration for  variations  in  ) ength-to-diameter  ratio(J./D  = 5,  7.5,  10,  12.5,  15) 
are  presented  in  APPENDIX  7.0. 
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Near  the  matched  pressure  point,  separation  should  not  be  a problem 
and  the  mass  flow  ratio  data  (triangles)  of  Figs.  3.3-1  to  3.3-4  should 
lie  on  the  pressure  data  (circles).  For  NL  = 1.75  and  2.50,  this  is 
indeed  true;  but  for  My  = 1.50  and  2.00,  there  is  a small  deviation. 
This  deviation  is  apparently  due  to  boundary  layer  growth  in  these  two 
nozzles  caused  by  the  constant-area,  exit  segments  which  were  added  to 
the  secondary  nozzle  contours  as  specified  in  Fig.  3.1-4. 

As  mentioned  previously,  the  one-dimensional  theory  of  Section  2.0 
employed  an  isentropic  recompression  model  for  prediction  of  the  upper 
limit  to  the  plane  of  supersonic- supersonic  operation.  In  Fig.  3.3-12, 
the  wall  pressure  distributions  near  the  upper  limit  point  in  Figs.  3.3- 
to  3.3-8  have  been  replotted  with  Pso  as  the  characteristic  pressure  in 
order  to  compress  the  scale.  The  initial  recompression  of  the  secondary 
stream  to  Mach  numbers  probably  less  than  1 (P(x)/P  < 0.5283)  as  shown 

in  Fig.  3.3-12  suggests  another  simple  model,  a "two  shock  model"  as 
illustrated  in  Fig.  3.3-13,  which  is  also  suitable  for  a one-dimensional 
analysis  and,  in  addition,  accounts  for  the  irreversibilities  in  two 
oblique  shocks.  Of  course,  this  model  does  not  allow  for  secondary  flow 
separation,  nor  would  it  apply  to  large  turning  angles  where  the  two 
oblique  shocks  are  replaced  by  a single,  normal  shock  wave. 

An  attempt  was  made  to  compare  the  one-dimensional  theory  of 
Section  2.C  with  the  extensive  sup  rsonic  wir.d  tunnel -injector  data  or 
Basel  and  Sinclair  [4].  Unfortunately,  these  investigators  did  not 
measure  the  secondary  stagnation  pressure  Pgo  or  include  the  static 
pressure  P at  the  injector,  which  makes  any  comparison  impossible. 
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Figure  3.1-2  Section  View  and  Specifications  of  the  Continuous-Slope 
Nozzle  and  Constant-Area  Mixing  Tube  Configurations 
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Figure  3.1-3  F.nlargeinent  of  a Typical  Continuous-Slope,  Supersonic- 
Supersonic  Nortzle 
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F.jector  Experiment  Flow  Diagram 
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Maximum  Compression  Characteristics  for  the  Constant -Area, 
Supersonic-Supersonic  Ejector  (M  = 1.75) 


Maximum  Compression  Characteristics  for  the  Constant-Area, 
Supersonic-Supersonic  Ejector  (M  = 2.001 
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Figure  3.3-5  Wall  Pressure  Distributions  for  the  Constant-Area, 
Supersonic-Supersonic  Ejector  at  Maximum 
Compression  Conditions  (M  = 1.50) 
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Figure  3.3-9  Exit  Mach  Number  Distributions  for  the  Constant 
Area,  Supersonic-Supersonic  F.jector  at  Maximum 
Compression  Conditions  (M  = 2.00) 
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(b)  With  Secondary  Separation 


Figure  3.3-10  Schematic  of  Secondary  Flow  Separation 
Induced  by  the  Primary  Flow 
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Figure  3.3-11  Variation  in  the  Secondary  Mach  Number  at  the 
Mixing  Tube  Entrance  with  Primary  Stagnation 
Pressure  and  Secondary  Nozzle  Reynolds  Number 
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Figure  3.3-12  Wall  Pressure  Distributions  for  the  Constant- 
Area,  Supersonic-Supersonic  Ejector  Near  the 
Upper  Limit  Point 
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Figure  3.3-13  Schematic  of  a "Two  Shock"  Model  for  the  Constant- 
Area,  Supersonic-Supersonic  Ejector 
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4.0  EJECTOR  OPTIMIZATION  AND  COMPARISON  OF  THE  CONSTANT- AREA, 
SUBSONIC-SUPERSONIC  AND  SUPERSONIC-SUPERSONIC  EJECTORS 

One  purpose  of  this  investigation  was  to  compare  the  performance  of 
the  constant-area,  supersonic-supersonic  ejector  with  that  of  the  constant- 
area,  subsonic-supersonic  ejector  as  applied  to  high-energy  chemical  laser 
systems.  Because  of  the  large  number  of  parameters  involved,  it  appears 
that  the  only  fair  comparison  of  the  two  ejectors  must  be  based  on  optimum 
data  for  a supersonic-supersonic  pumping  system  and  a subsonic-supersonic 
pumping  system.  To  this  end,  a method  was  developed  for  optimizing  these 
two  pumping  systems  within  a given  set  of  constraints.  This  method  is 
based  entirely  on  one-dimensional,  compressible  flow  theory  with  a limited 
use  of  empirical  data. 

4.1  RELATIONSHIP  OF  THE  CONSTANT-AREA  SUBSONIC-SUPERSONIC  EJECTOR 

TO  THE  CONSTANT-AREA  SUPERSONIC-SUPERSONIC  EJECTOR 

The  performance  characteristics  of  the  constant-area,  subsonic- 
supersonic  ejector  and  constant-area,  supersonic-supersonic  ejector  are 
quite  similar,  the  primary  difference  being  the  range  of  M . 

Given  yg  , yp  , Mwg /Mwp  , Tso/Tpo,  and  A^/A^,  the  performance  char- 
acteristics of  the  constant-area,  subsonic-supersonic  ejector,  as  shown 
schematically  in  Fig.  4.1-1,  are  described  by  a three-dimensional  surface 
with  axis  M^  , psi/pp0>  and  pm/ppo  . Since  the  secondary  entrance  forms 
a converging  nozzle,  Mgl  may  take  on  any  value  in  the  range  0 to  1 . For 
P > Psi  , the  primary  stream  expands  inside  the  mixing  tube  and  the 
secondary  stream  is  reaccelerated  to,  at  most,  sonic  conditions.  Should 
the  secondary  stream  be  thus  choked,  the  ejector  operation  becomes  inde- 
pendent of  P /P  as  indicated  by  two  faces  of  the  characteristic 

r M3  P 0 
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surface  which  lie  parallel  to  the  P^/Ppg  axis.  The  curve  .narking  the 

transition  from  the  P^a/Ppo  independent  to  dependent  regimes  is  known  as 

the  "break-off"  curve  and  is  considered  to  represent  the  optimum 

t 

conditions  for  ejector  operation. 

Given  yg  , Yp  . Mwg/Mwp,  Tgo/Tpo,  and  the  performance  char- 

acteristics of  the  constant-area,  supersoni c-supersonic  ejector,  as  shown 
schematically  in  Fig.  4.1-2,  are  also  described  by  a three-dimensional 
surface  with  axis  ML  , Pcl/PDn,  and  P,„/Pnn;  however,  the  ejector 

>5  1 i>  1 r 0 m3  P 0 

operation  is  restricted  to  only  a portion  of  the  surface.  So  long  as  the 
entering  secondary  stream  remains  supersonic,  as  assumed  in  Section  2.0, 
will  be  at  the  supersonic  design  value  as  produced  by  the  generating 
device,  e.g.  a laser  cavity,  and  the  ejector  is  confined  to  the  plane  of 
supersonic-supersonic  operation.  If  Pm/Ppo  is  increased  beyond  the 
limit  of  maximum  Pj^/PpQ  for  supersonic-supersonic  operation,  then  a 
normal  shock  wave  will  pass  into  the  secondary  entrance  and  M must 
undergo  a step  change  to  the  normal  shock  value  as  indicated  in  Fig. 
4.2-2(a).  If,  alternately,  Psl/Ppo  is  decreased  oeyond  the  upper  limit 
line  for  supersonic-supersonic  operation,  then  the  secondary  stream  is 
recompressed  in  the  mixing  tube  to  an  area  less  than  A*  for  the  super- 
sonic design  value  of  , and  must  drop  to  the  subsonic  design  value 
of  the  generating  device.  However,  if  Ps!/Pp0  is  now  increased,  the 
entering  secondary  flow  will  again  go  supersonic,  similar  to  a second- 
throat,  supersonic  diffuser,  with  Mgj  first  increasing  to  the  normal 
shock  value  before  a step  change  to  the  supersonic  design  value. 

4* 

'See  references  [1,8,10,11]  for  a complete  description  of  constant-area, 
subsonic-supersonic  ejector  operation. 
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The  similarity  in  performance  characteristics  of  the  constant-area, 
subsonic-supersonic  ejector  and  constant-ar~r , supersonic-supersonic 
ejector  becomes  even  more  apparent  as  the  supersonic  design  value  of  M 
in  Fig.  4 . 1 - 2 (a ) approaches  1 since  the  normal  shock  value  also  approaches 
1 and,  in  the  limit,  Fig.  4.1-2(a)  becomes  identical  to  Fig.  4.1-l(a). 

4.2  OPTIMIZATION  OF  THE  CONSTANT- AREA,  SUBSONIC-SUPERSONIC  AND 
SUPERSONIC-SUPERSONIC  EJECTORS 

Referring  to  the  high-energy,  chemical  laser  system  schematic  of 
Fig.  4.2-1,  the  laser  cavity,  stations  1 to  2,  has  no  influence  on  the 
pumping  system  other  than  establishing  the  flow  conditions  at  point  2. 
Hence,  the  laser  cavity  can  be  replaced  by  a supersonic  wind  tunnel  or  any 
other  device  producing  a uniform,  supersonic  stream.  A one-dimensional 
analysis  for  laser  cavity  flows  with  heat  addition  is  available  [1]  and 
was  included  in  a computerized  version  of  the  optimization  procedure. 

The  constant-area,  normal  shock  diffuser,  stations  2 to  3 of  Fig. 
4.2-1,  diffuses  the  entering  supersonic  stream  to  its  normal  shock  value. 
The  static  pressure  ratio  across  tho  duct  is  expressed  by 


VP2  ■ 


NSD 


V 


where  R^s D is  an  empirical  normal  shock  coefficient  and  ffy,  , M2 ) is  th 
usual  normal  shock  static  pressure  ratio  function.  Value:  of  T‘  D in.  tne 
range  0.75  <_  <_  1.25  are  commonly  used  for  parametric  studies  [ 1 j. 

The  static  pressure  rise  ac70ss  the  subsonic  diffuser,  stations  3 v_ 

4 of  Fig.  4.2-1,  is  given  by 


(P. 
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where  ri  is  an  empirical  diffuser  efficiency  and  (P  - P3)IDEAL  is  a 
function  of  Ys > . and  A4/A^.  As  a simplification,  n was  taken  from  a 

curve  fit  of  experimental  data  given  in  [17,18]  for  15°  conical,  subsonic 
diffusers  with 

H = 1.048992  - 0.027229  (A  /^) 

- 0.024461  (A4  /A3 ) 2 + 0.002685  (A4/A3)3  , 

1 • 0 i a4  / Aj  i 5 • 2 • 

This  empirical  relation  was  also  applied  to  the  subsonic  diffuser  at 
stations  7 to  8. 

Although  stagnation  pressure  losses  are  encountered  in  sudden 
enlargements,  these  losses  were  ignored  in  the  subsonic-supersonic  pump- 
ing system  of  Fig.  4.2-1. 

Finally,  the  one-dimensional  analysis  for  the  constant-area, 
supersonic-supersonic  ejector  of  Fig.  4.2-1  was  taken  from  Section  2.0, 
and  the  corresponding  analysis  for  the  constant-area,  subsonic-supersonic 
ejector  was  taken  from  Addy  and  Mikkelsen  [l]. 

Then  given  Y$ , M„ , R^,  A^,  yp , Mwp/Mw$,  T^/^,  A8/A?  , P&0/P2> 
and  P8/P2  for  the  subsonic- supersonic  pumping  system,  or  yg  , , yp , 

Mwp/Mws,  T60/T20  , A8/A7,  P6Q/P2,  and  Pg/P2  for  the  supersonic-supersonic 
pumping  system,  tho  optimum  is  considered  ."o  be  that  ejector,  as  specified 
by  M.  and  A„/A,  which  requires  the  minimum  WL/HL;  or,  alternately,  the 
minimum  P60/P2  when  Wp/Wg  is  specified.  In  general,  y , M2  , yp  , Mwp/Mwg, 
and  T60/T50  (or  T60/T20)  are  known  and  Pg/P2  and  P60/P2  (or  Wp/Wg)  are 
the  constraints  of  primary  interest. 

Now  consider  the  simplified  supersonic-supersonic  pumping  system  of 
Fig.  4.2-1  with  Ag/A?  - 1,  i.e.  with  no  subsonic  diffuser  and  y , yp , 
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MWd/MWcj  T,/T  and  H.  as  known  constants.  Then  given  M,  and  A, /A,  , 

upper  limit  point  values  of  P?/P2,  P60/P2  » and  Wp/Wg  are  obtained  from 
the  one-dimensional  theory  of  Section  2.0.  Figure  4.2-2  is  a plot  of 
M&  vs.  P?/P2  at  the  upper  limit  point  over  a range  of  A2/A&  for  this 
system.  If  the  compression  ratio  P?/P2  is  chosen  to  be,  say,  20,  then 
all  possible  combinations  of  and  A2/Afe  satisfying  this  constraint  are 
determined  by  the  verticle  dashed  line. 

Figure  4.2-3  is  a cross  plot  of  Pfc0/P2  and  Wp/Wg  vs.  A2/Afe  for  the 
combinations  of  M&  and  A2/A&  of  Fig.  4.2-2  satisfying  the  constraint  on 
P_/P* . From  this  graph  it  is  apparent  that  P /P  varies  inversely  with 
Wp/Wg  and  that  if  an  upper  limit  is  established  for  P60/P2>  say,  at  1000, 
then  the  minimum  value  of  W /W  , in  this  case  W /W  = 3.2,  will  occur  at 
that  upper  limit.  The  optimum  ejector  is,  therefore,  fully  determined 
since  = 4.4  for  P60/P2  = 1000  and  Mfc  = 4.4  from  Fig.  4.2-2.  If, 

for  example,  an  upper  limit  of  6 is  set  for  W /W  , then  the  minimum  value 

P o 

of  P60/P2»  in  this  example  P60/P2  = 217,  will  occur  at  \/\  = 1,  corre- 
sponding to  the  upper  limit  value  of  Wp/Wg  in  Fig.  4.2-3  and  = 3.7 
from  Fig.  4.2-2. 

This  procedure  for  optimizing  the  supersonic-supersonic  pumping 
system  of  Fig.  4.2-1  is  not  significantly  altered  by  the  addition  of  a 
subsonic  diffuser  at  stations  7 to  8.  Figures  4.2-2  and  4.2-3  are  still 
sufficient;  however,  the  overall  compression  ratio  P./P,  is  substituted 
for  P./P,  . 

The  optimization  procedure  for  the  subsonic-supersonic  pumping  system 
of  Fig.  4.2-1  is  nearly  identical  to  the  procedure  outlined  above  for  the 
supersonic-supersonic  pumping  system,  with  Figs.  4-2-4  and  4.2-5 
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corresponding  to  Figs.  4.2-2  and  4.2-3,  respectively;  however,  one  addi- 
tional variable  is  involved  since  Mg  may  take  on  any  value  in  the  range 
0 <_  Mg  <_  1.  The  result  is  that  Figs.  4.2-4  and  4.2-5  must  be  reproduced 
for  all  values  of  Mg  which  give  solutions  satisfying  the  constraints  or. 
P7/'P2  anc*  P60/^2  (°r  alternately  Wp/W  ).  Once  this  is  accomplished,  the 
minimum  values  of  Wp/Wg  (or  alternately  P60/P2)  can  136  cross  plotted  vs. 

M.  to  find  the  absolute  minimum  as  shown  in  Figs.  4.2-6  and  4.2-7.  This 
process  completes  the  ejector  specification  since,  for  example,  W /W 
and  ^ from  Fig.  4.2-6  determine  Ag/A&  from  Fig.  4.2-5,  which  in  turn 
specifies  from  Fig.  4.2-4. 

As  an  example  cf  the  optimization  procedure  for  a subsonic-supersonic 
pumping  system  of  Fig.  4.2-1  with  1^  , A /A,  and  A /A?  equal  to  1,  let 
Mg  - 0.9  as  in  Figs.  4.2-4  and  ‘.2-5  with  constraints  of  P?/P2  = 20  and 
P60/P2  - 1000.  Since  P60/P2  varies  inversely  with  Wp/Ws  in  Fig.  4.2-5, 
the  minimum  value  of  W /W_  is  2.8  at  Ac /A,  - 3.8  and  M^  =4.5  from 
Fig.  4.2-4.  Repeating  this  process  for  all  suitable  values  of  Mc  and 
plotting  the  results,  shown  in  Fig.  4.2-6,  it  can  be  seen  that  the  abso- 
lute minimum  value  of  Wp/Wg  occurs  at  Mt  = 0.9  and  the  optimization  is 
complete. 

If  in  the  preceding  example  an  upper  limit  of  hp/Wg  = 6 was  sel^ted 
rather  than  the  constraint  on  P../P. , then  for  M,  = 0.9,  the  minimum  value 
of  P,„/P„  would  be  105  at  Ac/A,  = 0.8  from  Fig.  4.2-5  and  M,  = 2.8  from 
i ig.  4.2-4.  Now  repeating  this  process  for  all  values  of  Mg  yielding 
solutions  which  satisfy  the  system  constraints  and  plotting  the  resuits 
shown  in  Fig.  4.2-7,  it  can  be  seen  that  the  absolve  minimum  value  of 
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P,./P„  occurs  at  M,  =0.7,  rather  than  at  = 0.9,  and  that  the  solution 
at  Mg  = 0.9  is  merely  one  step  in  the  optimization  process. 

While  this  graphical  method  of  ejector  optimization  is  necessarily 
quite  tedious,  the  entire  process  is  readily  adaptable  to  standard  com- 
puter programing  procedures.  Some  care,  however,  must  be  taken  in  the 
selection  of  Mg,  Ag/Ag,  and  (or  /A  and  M& ) since  only  certain  com- 
binations will  yield  solutions  satisfying  any  one  system  constraint. 

4.3  COMPARISON  OF  OPTIMUM  CONSTANT-AREA,  SUBSONIC-SUPERSONIC  AND 

SUPERSONIC-SUPERSONIC  EJECTOR  DATA 

Three  sets  of  optimum  data,  as  listed  in  Tables  4.3-1,  4.3-2,  and 
4.3-3,  were  calculated  by  the  procedure  of  Section  4.2  for  comparison  of 
a supersonic- Jupersonic  pumping  system  with  a subsonic-supersonic  pumping 
system  as  shown  in  Fig.  4.2-1.  The  first  two  cases  represent  typical 
high-energy,  chemical  laser  system  data  while  the  third  case  demonstrates 
a supersonic  wind  tunnel  application.  In  each  case,  the  appropriate 
variable,  either  Wp/Wt,  or  P& Q / ?2  , was  minimized  for  the  supersonic- 
supersonic  pumping  system  and  for  the  subsonic-supersonic  pumping  system 
with  = 1.0,  0.85,  and  0.75.  The  exit-to-entrance  area  ratio  for  all 

the  subsonic  diffusers  was  arbitrarily  set  at  2.0.  The  optimum  ejecto. 
for  each  case  is  specified  by  M and  A^/A  The  area  ratio  Ag/A2  was 
also  tabulated  as  an  indication  of  the  overall  size  of  the  pumping 
systems . 

In  selecting  data  for  the  hypothetical,  supersonic  wind  tunnel  con- 
figuration of  Table  4.3-3,  it  was  assumed  that  the  secondary,  or  test 
section,  nozzle  and  the  ejector  primary  nozzle  would  both  be  supplied 
from  a common  source  at  the  same  stagnation  pressure,  in  this  case 
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790,8  kPa,  and  that  the  ejector  system  would  pump  to  atmospheric  condi- 
tions. This  assumption  precludes  a minimization  of  P&0/P2  since  the 
compression  ratio  Pg/P2  would  be  unknown. 

The  optimization  method  of  Section  4.2  used  only  theoretical, 
supersonic- supersonic  ejector  data  computed  at  upper  limit  point  condi- 
tions. In  defense  of  this  procedure,  Cases  No.  1 and  2 were  also 
calculated  for  a supersonic- supersonic  pumping  system  at  matched  pressure 
conditions.  A matched  pressure  calculation  was  unrealistic  for  the 
supersonic  wind  tunnel.  Case  No.  3,  since  this  condition  would  only  be 
satisfied  if  the  secondary,  or  test  section,  nozzle  and  primary  ejector 
nozzle  (operating  at  equal  stagnation  pressures)  had  identical  design 
Mach  numbers . 

Examination  of  the  optimum  data  for  Cases  No.  1,  2,  and  3 leads  to 
several  important  conclusions: 

1.  For  a fixed  value  of  P,  /P„,  the  minimum  value  of  W/W„  for  the 
supersonic-supersonic  pumping  system  lies  between  the  minimum  value  of 
.Vp/Ws  for  the  subsonic-supersonic  pumping  systems  with  =1.0  and 

0.85.  Since  a value  of  0 85  is  most  realistic,  it  appears  that 

the  supersonic-supersonic  pumping  system  gives  the  best  performance. 

For  a fixed  value  of  Wp/Wg,  the  minimum  value  of  P60/P2  for  the 
supersonic-supersonic  pumping  system  is  equal  to  or  greater  than  the  mini- 
mum value  of  P60/P2  for  the  subsonic-supersonic  pumping  system  with 
Rjj  - 0.75  indicating  that  the  performance  of  the  subsonic-supersonic 
pumping  system  is  superior. 

Clearly,  the  only  conclusion  to  be  drawn  from  this  limited  data  is 
that  the  selection  of  one  pumping  system  over  the  other  on  the  bases  of 
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mass  flow  ratio  and  primary  stagnation  pressure  depends  entirely  on  the 
system  constraints  and  that  no  broad  statement  can  be  made  as  to  one 
pumping  system  performing  better  than  the  other. 

2.  Upper  limit  point  conditions  are  preferred  over  matched  pressure 
conditions  in  the  optimization  of  supersonic-supersonic  pumping  systems 
since  the  resultant  system  requires  much  lower  values  ot  Wp/W  and 
P60/P2.  The  same  conclusion,  though  not  proven,  should  apply  to  all  values 
of  P&/P2  between  the  matched  pressure  and  upper  limit  points  since  P?/P2 
varies  line' rly  with  W /IV  and  P,  /P„  (see  Section  2.0). 

3.  Based  on  the  area  ratio  A„/A„,  that  pumping  system  with  the  best 
performance,  or  the  smallest  !Vp /Wg  or  P60/P2  requirements,  is  also 
physically  the  smallest  system. 

4.  The  normal  shock  coefficient  has  a significant  influence  on  the 
performance  of  subsonic-supersonic  pumping  systems  and  points  to  super- 
sonic diffuser  design  as  an  area  for  careful  attention. 

5.  The  comparison  of  a constant-area,  subsonic-supersonic  ejector 
with  a constant-area  supersonic-supersonic  ejector  on  the  bases  of 
primary  Mach  number,  Mfe , and  mixing  tube  area  ratio,  A?/A6>  alone  has 
little  if  any  value. 

4.4  COMPUTER  PROGRAMS 

Two  computer  programs  were  developed  to  conduct  the  foregoing  opti- 
mization and  system  studies.  The  optimization  program,  CLGDOP,  and  the 
systems  program,  CLGPSP,  along  with  sample  input  and  output  data  are 
presented  in  detail  in  APPENDICES  7.4  and  7.5,  respectively. 


independent 


Pms/Ppo  ' 


Break-oiff 

curve 


independent /N  D /B 
Js  / | "nij/Hpo 

^ / \ dependent 


Pmj/Ppo 


(a)  Ejector  Characteristic  Surface 


PO. Ipo 


Primary  MPI,Pp,,Ap, 


Secondary  | Msl,P8l,Asl 
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Figure  4.1-1  Constant-Area,  Subsonic-Supersonic  Ejector 
Performance  Characteristics 
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(a)  Ejector  Characteristic  Surface 


Pro  «^?0 


f b)  Ejectoi  Configuration  and  Notation 


Figure  4.1-2  Ccnstant-Area,  Supersonic-Supersonic  Ejector 
Performance  Characteristics 
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Figure  4.2-4  Typical  Break-Off  Data  for  Optimization 
of  a Subsonic-Supersonic  Pumping  System 
(M6  vs.  P?/P2) 


I 


Figure  4 . 2 
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-5  Typical  Break-Off  Data  for  Optimization 
of  a Subsonic-Supersonic  Pumping  System 
(W„/K'  and  P,„/P„  vs.  A/'AJ 


85 


4 


Table  4.3-1  Optimum  Chemical  Laser  System  Data,  Case  No.  1 

(a)  Minimum  W /W 
p'  s 


CASE  NO.  1(a) 

Optimum  Chemical  Laser  System  Data 

for  Minimum  W /W 
p s 


= 1 562 

Mwp  /Mwg 

= 1.684 

^p 

= 1,340 

K /K 

8 7 

= 2.000 

M2 

= 2.180 

P /P 
60 ' 2 

= 1937 

P./P, 
8 2 

= 28.06 

Supersonic-Supersonic  Pumping  System 
at  Upper  Limit  Point  Conditions 

WT20  ’ °'761 

W /V 

P S 

M6 

V\ 

am 
8 2 

4.408 

4.690 

5.028 

2.496 

Subsonic-Supersonic  Pumping  System 
at  Break-Off  Conditions 


T /T 

*6  0'  5C 

= 0.761  A4/Aj 

= 2.000 

"nsd 

w /IV 
P s 

M6 

V\ 

A. /A. 
8 2 

1.00 

3.747 

4.734 

4.828 

2.125 

0.85 

4.792 

4.858 

4.136 

2.614 

0.75 


5.777 


4.955 


3.695 


3.079 
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Table  4.3-1  (Cont.) 

(b)  Minimum  P60/P2 

CASE  NO.  1(b) 

Optimum  Chemical  Laser  System  Data 

for  Minimum  P, , /P„ 

6 0 2 

Y_  = 1.562  Mw  /Mw  = 1.684 

s PS 

Yn  = 1.240  Aa /A„  = 2.000 

'p  8 7 

M2  = 2.180  Wp/Ws  = 10.00 

Pg/P2  = 28.06 

Supersonic-Supersonic  Pumping  System 
at  Upper  Limit  Point  Conditions 

T60/T20  * °'761 

P60/P2  M6  VA6  A8/A2 

318.6  3.797  1.705  4.838 


Subsonic-Supersonic  Pumping  System 
at  Break-Off  Conditions 


1 6 0 7 ^5  0 

= 0.761  A 4/Aj 

= 2.000 

^USD 

?60/P2 

M6 

v\ 

VAa 

1.00 

110.2 

2.660 

1.729 

4.319 

0.85 

174.0 

3.047 

1.903 

4.515 

0.75 

270.6 

3.410 

2.085 

4.665 
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Table  4.3-1  (Cont.) 

(c)  Minimum  Wp/Ws  at  Matched  Pressure  Conditions 


CASE  NO.  1(c) 


Optimum  Chemical  Laser  System  Data 

for  Minimum  W /W„ 
p'  s 


= 1.562 

Mw/Mw  = 1.684 
p ' s 

Yp 

= 1.340 

= 2.000 

M2 

= 2.180 

P60/P2  " 1937 

P8/P2  = 28.06 


Supersonic-Supersonic  Pumping  System 
it  Matched  Pressure  Conditions  (P, /P_  = 1.0) 

Q l 

\J\*  - 0-761 

V\  «6  V\  VA2 


i 

f 


9.665 


5.853 


1.654 


5.058 
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Table  4.3-1  (Cont.) 

(d)  Minimum  P,  /P„  at  Matched  Pressure  Conditions 

o 0 2 


CASE  NO.  1(d) 

Optimum  Chemical  Laser  System  Data 

for  Minimum  P,  / P„ 

60  2 


Ys 

= 1. 

.562 

Mw  /Mw„ 
P s 

= 1.684 

Yp 

= 1. 

.340 

A0/An 

8 7 

= 2.000 

m2 

= 2. 

.180 

w /wo 

P s 

= 10.00 

P„/P„ 
8 2 

= 28.06 

Supersonic-Supersonic  Pumping  System 
at  Matched  Pressure  Conditions  (P6/p2  = 1.0) 

VLo  ■ t>-761 

P60/P2  H6  V\,  \'\ 

1836  5.807  1.623  5.209 
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Table  4.3-2  Optimum  Chemical  Laser  System  Data,  Case  No.  2 

(a)  Minimum  Wp/Ws 


CASE  NO.  2(a)  j 

Optimum  Chemical  Laser  System  Data  I 

for  Minimum  NL/W  i 

p s I 


Ys 

= 1.562 

Mw„  /Mw 
p s 

•*  1.478 

Yr 

= i.340 

VA-7 

- 2.000 

M2 

= 2.230 

P60/P2 

= 1430 

P„/P„ 
a ' 2 

=»  20.71 

Supersonic-Supersonic  Pumping  System 
at  Upper  Limit  Point  Conditions 

V5o  * e.817 

\ V\  VA2 

2.228  4.055  1.053  2.209 


Subsonic-Supersonic  Pumping  System 
at  Break-Of?  Conditions 


T /T  = 
to'  so 

0.807  A,  /A,. 

= 2.000 

p 

NS  D 

W D/Wc 
P s 

M6 

VA6 

v-s 

1.00 

1.968 

4.476 

7 . 000 

1.952 

0.85 

2.559 

4.596 

5.838 

2.375 

0.75 


3.108 


4.690 


5.128 


2.771 
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Table  4.3-2  (Cont.) 
(bj  Minimum  P,  /?„ 

> vi  0 2 


CASE  NO.  2(b) 

Optimum  Chemical  Laser  System  Data 


for 

Minimum  P,  /P„ 
60  2 

Ys  * 1-562 

Mw  /Mw  = 1.478 
p s 

Yp  = 1-340 

Ag /A?  = 2 . 000 

M2  = 2.230 

Wp/Wg  = 6.000 

P8/P2  = 20.71 

Supersonic-Supersonic  Pumping  System 
at  Upper  Limit  Point  Conditions 

■VT20  = °-8»7 

V*C  VA2 

123.0  3.127  1.801  4.497 


Subsonic-Supersonic  Pumping  System 
at  Break-Off  Conditions 


or3 

O 

cT3 

O 

= 0.807  A4/A3 

= 2.000 

^30 

P /P 

M6 

VA6 

VA2 

1.00 

56.38 

2.223 

■*T 

r? 

CO 

r— 1 

3.951 

0.85 

84.21 

2.568 

2.C47 

4.160 

0.75 

12S.7 

2.894 

2.262 

4.323 
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Table  4.3-2  (Cont.) 

(cj  Minimum  bp/W$  at  Matched  Pressure  Conditions 


CASE  NO.  2[c) 

Optimum  Chemical  Laser  System  Data 

for  Minimum  W /W„ 
p s 


<s 

= 1.562 

Mv;p/Mws  = 1.47b 

Yp 

= 1.340 

= 2.000 

M2 

= 2.230 

P60/P2  ■=  1430 

Pg/P?  = 20.71 


Supersonic-Supersonic  Pumping  System 


at  Matched 

Pressure  Conditions  (P, /P„ 

o 2 

T60/T20  = °-?07 

= 1.0) 

Vws 

M6 

VS 

\/A2 

3.920 

5.593 

2.302 

3.536 

Table  4.3-2  (Cont.) 

(d)  Minimum  P /P  at  Matched  Pressure  Conditions 


CASE  NO.  2(d) 

Optimum  Chemical  Laser  System  Data 
for  Minimum  Pfi0/P2 


Ys  - 1-562 

Mw  /Mw  = 1.478 
p s 

Yp  = 1-340 

O 

o 

o 

(N 

II 

<r 

oo 

< 

M2  = 2.230 

Wp/Ws  = 6.000 

vp2  - 

20.71 

Supersonic-Supersonic  Pumping  System 
at  Matched  Pressure  Conditions  (P, /P_  = 1.0) 

D 2. 

o 

o 

= 0.807 

P /P  M 

60 ' 2 6 

A8/A2 

633.0  4.934 

1.677  4.955 
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Table  4.3-3  Optimum  Supersonic  Wind  Tunnel  Data 


CASE  NO.  3 

Optimum  Supersonic  Wind  Tunnel  Data 


for  Minimum  W?/Wg 

Ys 

= 1, 

.400 

Mw„  /Mw 
p s 

= 1.000 

= 1, 

.400 

VA 

= 2.000 

= 5, 

.000 

P60/P2 

= 529.1 

Co 

*0 

to 

)) 

O 

-4 

OO 
►— * 

Supersonic-Supersonic  Pumping  System 
at  Upper  Limit  Point  Conditions 


o 

O 

= 1.000 

W /Wc 
P s 

M, 

O 

V\ 

AM 
& 2 

1.810 

2.063 

8 . 758 

2.258 

Subsonic-Supersonic  Pumping  System 
at  Brea!:-0ff  Conditions 


T /T 

60'  50 

= 1 000 

= 2.000 

^NSD 

w /w„ 
P ' s 

M6 

VA2 

1.00 

1.349 

2.811 

4.450 

1.700 

0.85 

1.852 

2.918 

3.669 

2.130 

0.75 

2.336 

5.001 

3.217 

2.545 
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5.0  CONCLUSIONS 

As  a result  of  this  preliminary  theoretical  and  experimental  analysis 
of  the  constant-area,  supersonic-supersonic  ejector,  the  following  con- 
clusions may  be  drawn. 

1.  A one-dimensional  theory  was  developed  which  predicts  the 
performance  characteristics  of  all  constant-area,  supersonic-supersonic 
ejector  configurations. 

2.  Due  to  the  simplified  analysis  of  the  constant-area,  mixing 
section  and  inviscid  interaction  region,  the  present  one-dimensional 
theory  is  particularly  well-suited  to  parametric  evaluations  of  constant- 
area,  supersonic-supersonic  ejector  performance. 

3.  The  present  one-dimensional  theory  predicts  maximum  compression 
ratios  which  are  15  to  21  percent  greater  than  experimental  measurements. 

4.  The  constant-area,  supersonic-supersonic  ejector  is  particularly 
susceptible  tc  secondary  flow  separation  which  requires  a more  sophisti- 
cated method  of  analysis  than  the  present  one-dimensional  theory. 

5.  Continued  experimentation  is  needed  to: 

a)  Establish  a sufficient  length  for  constant-area,  supersonic- 
supersonic  ejector  mixing  tubes. 

b)  Obtain  supersonic-supersonic  ejector  data  over  a wide  range 
of  operation,  particularly  under  simulated  high-energy, 
chemical  laser  system  conditions. 

c)  Verify  the  present  one-dimensional  theory  for  plane,  two- 
dimensional  ejector  configurations  which  are  more  consistent 
with  current  high-energy,  chemical  laser  designs. 
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d)  Develop  more  advanced  theoretical  models,  including  a "two- 
shock"  model,  through  flow  visualization  studies  of  the 
mixing  and  interaction  phenomena. 

e)  Evaluate  the  effects  of  variable  mixing  tube  wall  profiles 
which  have  proved  advantageous  in  subsonic-supersonic 
ejector  development. 

6.  The  constant-area,  supersonic-supersonic  ejector  is  an  extension 
of  the  constant -area,  subsonic- supersonic  ejector  with  similar  character- 
istic surfaces. 

7.  The  constant-area,  subsonic-supersonic  and  supersonic-supersonic 
ejectors  must  be  compared  on  an  optimum,  overall  pumping  system  basis. 

8.  Present  one-dimensional  theories  for  both  the  constant-area, 
subsonic-supersonic  and  supersonic-supersonic  ejectors  may  be  incorpor- 
ated in1'  a single  high-energy,  chemical  laser  system  optimization 
procedure. 

9.  Conclusions  as  to  the  relative  performance  of  optimum  subsonic- 
supersonic  and  supersonic-supersonic  pumping  systems  are  included  in 
Section  4.3;  but  in  summary,  a supersonic-supersonic  pumping  system  has 
the  potential  for  improved  performance  over  that  of  a subsonic-supersonic 
pumping  system  on  the  basis  of  mass  flow  ratio  and  primary  stagnation 
pressure;  however,  this  depends  on 

a)  The  high-energy  chemical  laser  system  constraints, 

b)  The  source  of  disagreement  between  the  present  one-dimensional , 
supersonic-supersonic  ejector  theury  and  experiment  as  noted 

in  item  3,  preceding,  and 


96 

c)  An  adequate  solution  to  the  secondary  flow  separation  problems  en- 
countered during  high  compression  ratio  operation  of  the  constant- 
area,  supersonic-supersonic  ejector. 
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7.2  CONSTANT- AREA  SUPERSONIC-SUPERSONIC  EJECTOR  COMPUTER  PROGRAM 
7.2.3  Computer  Program (CASSE) 
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7.2.1  CASSE  (Cont.) 
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7.2.1  CASSE  (Cont.) 
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7.2.1  CASSE  (Cont.) 
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7.2.2  CASSE  Sample  Inp^it 


I NPUT  GS,  GP,  MWSMWP,  TS OTPO , MS  1 , MP 1 , AS  1 AP 1 
1 . 562,  1.34,  0.593912,  1.31341,  2.18,  4.690,  4.02823 

INPUT  RD, ERROR, NDATA 
0.75,  5.0E-06,  10 


CONSTANT-AREA  SUPER  SON  I C — SUPERS ON I C EJECTOR 


7.2.3.  CASSE  Sample  Output 
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7.2.3  CASSE  Sample  Output  (Cont 
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7.2.3  CASSE  Sample  Output  (Cont.) 
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7.2.3  CASSE  Sample  Output  (Cont.) 
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7.3  CONSTANT-AREA  SUPERSONIC-SUPERSONIC  EJECTOR  PARAMETERS 
COMPUTER  PROGRAM 

7.3.1  Computer  Program  (CASSEP) 
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7.3.1  CASSF.P  (Cont.) 
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INPUT  INDEPENDENT  VARIABLE  - TYPE  "H"  FOR  HELP 
H 


SELECT  THE  INDEPENDENT  VARIABLE  FROM  THE  FOLLOWING 
CATALOG  OF  EJECTOR  PARAMETERS: 

GS  /GP  , MWSMWP , TS  OTPO , MS  1 ,MP1  ,AS1AP1,RD 

INPUT  INDEPENDENT  VARIABLE  - TYPE  "H"  FOR  HELP 
AS 1 API 

INPUT  GS  ,GP  ,MWSMWP, TS0TP0,MS1  ,MP1  , RD 
ERROR 

1.562,  1.34,  0.593912,  1.31341,  2.18,  4.690,  0.75,  5.0E-06 

INPUT  LOW  VALUE,  INCREMENT,  AND  HIGH  VALUE  OF  AS1AP1 

1.0,  1.0,  10.0 
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7.3.3 


CASSEP  Sample  Output  (Cont.) 
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7.3.3  CASSlcP  Sample  Output  (Cont.) 
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7.3.3 


CASSEP  Sample  Output  (Cont.) 
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7.3.3  CASSEP  Sample  Output  (Cont.) 
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7.3.5  CASSEP  Sample  Output  (Cont.) 
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7.4  CHEMICAL  LASER  GAS  DYNAMICS  OPTIMIZATION  COMPUTER  PROGRAM 
7.4.1  Computer  Program 
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7.4.1  CI.GDOP  (Cont.) 
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7.4.1  CGLDOP  (Cont.) 
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7.4.2  CLGDOP  Sample  Case  No.  1 


CLG'iOP  SAMPLE  INPUT  DATA:  URGE  NO.  1 

INPUT  DATA  P UP  LASER  UR|IITV  FiHALYSlS  BY  NAMELIST. 

CUPPEHT  "ALOES  RPE: 

; C A* 1 

.J : - 1 . 5680 O < Ml  =5.  OO0O0 1 RPR  1 =8 . 55784 < C 0 1 = 1 . OOOOO  - 
C OP-  O . O000O  ? 1 1PT S=8  1 s 7 
- C Pi* 1 N 1=8.83*  RPRl-1 . Oj  C 0 1 — 0 • O j NFTS=57 

INPUT  THE  "API  ABLE  TO  BE  MINIMIZED  FROM  THE  FOLLOWING  LIST 

■P60Pi"  DIMENSIONLESS  PPIMRPY  STAGNATION  PPESSUPE 
1 NPWS"  PP I MAF'Y-TO-SEC OMDAPY  MASS  FLOW  RATIO 


! *nnpi 

1 1 iFUT  SYSTEM  CONSTPAINTS  BY  NAMELIST. 

CUPPEHT  i 1 ALLIES  ARE: 

i c or  1st  s 

F 8P 1 =57 . 1 634  J WP WS= l . 0O0O0  j 7 
i C0MST8  PSP 1 =£0 . 7141 j WF'W8=6.07 

INPUT  THE  EJECTOR  MODEL  FROM  THE  FOLLOWING  LIST: 

i.hE'1  CONSTRNT-APER  EJECTOR 

SSE " SUPERSONIC -SUPERSON I C EJECTOR 

C SE 

INPUT  DATA  FOR  EJECTOR  ANALYSIS  BY  NAMELIST. 

CURRENT  "ALOES  ARE: 

E JECT8 

GF'= 1 . 34O00 » MWPMWS= 1 . 68 375 .«  T 60T 80=0 . 76 1 376 r A8A7=8 . OOOOO 
7E JEC T8  MWPMWS= 1 . 478 1 5 * T60T80=0 . 8073 1 37 
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7.4.2  CGLDOP  Sample  Case  No.  1 (Cont.) 
CLGBOP  SAMPLE  OUTPUT  BATA:  CASE  MO. 1 


SEARCH  IiRTR: 


NTT  PE 

MSI 

NPMS 

DNS  1 

!•  H I L 

1 

0. 1O0000E+O1 

0. 5SS735E+01 

0 . 1 00000E-*  00 

NO 

CL 

0. S0000OE+00 

O. 5776S4E+81 

0. 100000E+0O 

NO 

L 

0 . 3O0000E+00 

0. 577978E+01 

0. 100000E+00 

NO 

- 

0. 80O000E+00 

0 . 5776S4E-+  0 1 

0. 100000E+00 

NO 

HIGH  EMEPGV  CHEMICAL  LASER  SYSTEM  SIMULATION 
OME-D I MENS  I OMfiL  ANALYS I S 

A.L.  ADDY 
C.D.  Mil  I ELSEM 
H.P.  SANDBERG 

1 JANUARY  76 

MECHANICAL  ENGINEERING  DEPARTMENT 
UN I "EPS I TY  OF  ILLINOIS  AT  UPBANA-CHAMPAIGN 
IJPBANAj  ILLINOIS  61  SOI 

CAE  SOLUTION  FOP  MINIMUM  WPWS 

SYSTEM  DATA: 

POINT  1 LASER  CANITY  ENTRANCE  CONDITIONS 

Ml  - 0.500000E+01  GS  0.  I562OOE+01 

PIP 10  - 0.30635SE-0S  TIT 10  - O. 12461 1E+00 

POINT  2 LASER  CAI,ITY  EMIT  AND  NORMAL  SHOCT  DIFFUSER 
ENTRANCE  CONDITIONS 

GS  - 0. 1 562OOE+0 1 

MS  - 0 . S 1 SGOOE+0 1 A2A1  = 0.255724E+01 

PSPl  = 5! 2036S2E+O1  T£T l - 0. 515579E+01 

P3OP10  = 0. 65864  IE-01  TS0T10  = 8. 15000OE+O1 


s® 
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7.4.2  CGLDOP  Sample  Case  No.  J (Cont.) 


POINT  3 ? IOF;MhL  'SHOD  DIFFUSE?  EXIT  RNI  -SUE SONIC 

DIFFUSE?  ENTRANCE  CONDITIONS 

GS  = 0.  3 SSSOOE+O 1 

H3  = 0 . 57 1 SSSE+88  FlSfil  - 0.  £557c'4E+01 

F Jp  1 — li.851?c4E+Cil  T 3T 1 — O.  1 lUcTTE+Uc 

P EOF1 10  - 8 . SScSTSE -8 1 T SOT  18  = 3.  158888E+01 

POINT  4 SUBSONIC  E I F FUSEE  EXIT  RHD  SUDDEN  fHLfiPGEMENT 
ErrPftUCE  CONDITIONS 

GS  = 8.  1 jCcOuE  *-81 

N4  - 0.  c'4S"?4 lC+08  R4Rl  - 8.511448E+81 

F4P1  = 0.  10£353E+0c  T4T1  - 8.  1 1S44SL+8S 

F48P 18  - 8 . 3£7'?S?E-8  1 T48T16  - 0.  150000E+01 

KOINI  5 CCMSTflMT-fiRtfi  EJECT  OF  SEC.'MDRRY  NOZZLE  EXIT 
c or  id  n ions 

GS  - 0. 1 SSS08E+8 J 

N5  = 8,  SOOO8OE+80  Rf.fi  1 = 0 . SS7 1 OSE-t  0 1 

PfF'l  - 0 . SOf SSSE+0 .1  TfTl  * 0.  98884 5E +81 

F'fOPlCi  = 0.  3S79b9E-6i  1 SOT  18  - 8.  158888E+81 

POINT  8 CONSTRNT-RPER  EJECTOR  PPIMftPV  NOZZLE  EXIT 
CONDITIONS 

GP  - 0.  1 S4OO0E+O1  NUF'MNS  ~ 0.  lbS375E-*01 

NF'NS  - 0.  577*84 E+81 

lit.  - 0.495558E+81  fibfi i --  0.  1 Obf 46E-+ 01 

F'bOPl  =•  0.  394478E+84  TOOT  3 - 0 91b789E+81 

F'SOPIO  = 0.  lc'0751E+8c  T68T10  - 0.  1 HSObE+OJ 

POINT  7 C ONSTRHT- flREfi  EJECTOR  EXIT  RHD  SfJB SONIC 
DIFFUSER  ENTPRNCE  CONDITION:. 

CM  - 0. 1 37388E+81  MMMMWS  - 0. X58944E+01 
MMWS  — O . S 77 Ss4E+8 1 

M7  - 8 . 4 355Z8E+08  R7R 1 - 0 . 3937  3 3E+8 1 

P7F'l  - 0.  5i3050bE+8£  T7T 1 = 0. 93S738E-+ 01 

P78P16  = 8. 3S18S7E+80  T70T18  = 0. 1 30309E+01 

POINT  S SUBSONIC  DIFFUSER  EXIT  CONDITIONS 


CM 

MI'IWS 

- O. 3 373S8E+81 
= 8. 6776S4E+81 

MWMMWS  = 

U.  1 f c‘944E+U  1 

MS 

= 0. 1 99607E+00 

RSFll 

Li.  7o74bbE+0 1 

PSP  1 

= 0. 571 155E+02 

T ST  1 

O.  95S6S9E+01 

POOR 1 0 

~ 8. 179667E+00 

TS8T18  = 

y.  1 cUSUSE+U 1 
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7.4.2  CGLDOP  Sample  Case  No.  1 (Cont.) 


LASER  CAM I TV  DATA: 


ARTS 

= R1 

CPI 

= O.1O000OE+01 

COR 

= 

0. 00OO0OE+00 

GS 

= O. 1 5SR00E+0 1 

ARA1 

s 

0. R557R4E+01 

Ml 

- 0 . 5O0000E+0 i 

MR 

— 

0.R1S008E+O1 

PRF‘1 

= 0 . RO3SSRE+0 1 

TRT1 

= 

0. 515579E+01 

PRO!-'  1 0 

= O. S58S41E-81 

TR0T10 

=: 

0. i J0Q00E+0 1 

HOPMAL 

SHOCK  DIFFUSER 

DATA: 

GS 

- O. 15SRO0E+01 

F:NSD 

= 

O. 750C 30E+0O 

MR 

- 0. R1S00OE+01 

M3 

= 

0 • 57 1 Sc'yE+00 

PSP  R 

— U.41 S 1 S4E +0 1 

T3T2 

=• 

0. R13SS9E+01 

PS0PR0 

= 0.505397E+60 

T30TRO 

r: 

0. 100000E+01 

SUBSONIC  DIFFUSER  DATA: 

FSD 

= 0 . 9S5RS 1 E+00 

GS 

- 0. 1 5S2OOE+0 1 

A4A3 

0 . R00000E+O ] 

MS 

- O . 57 1 SRSE+0O 

M4 

— 

0 . R4R94 1 E+00 

P4P3 

- 0. 1 £0 1 7 RE+0 1 

T4T3 

= 

0 . 1 U7 407E+O 1 

P40P30 

= O. 985RS1E+O0 

T40T38 

0. 100000E+01 

CO! ISTAMT-AREA  EJEC TOP 

DATA: 

GS 

- 0. 1 SSRO0E+01 

GP 

= 

0. 1340OOE+01 

GM 

= 0. 1 373R8E+0 1 

MWPMWS 

= 

O.  1 S8375E-+  0 1 

A7AS 

= O. 3S9544E+01 

HPWS 

= 

0 . 577 SS4E+0 1 

M6 

= O . 495558E+0 1 

M7 

= 

O . 43 55 ROE  ■+ 00 

F'SOF'SO 

— U . OSS 1 73E+03 

T S0T50 

0 . 7S 1 37 6E+0y 

P7P50 

•=  O. 485S14E+01 

T7T50 

O. 774S3SE+00 

F7OP50 

= 0. 55R147E+01 

T70T50 

O . SOROS 1 E+00 

1 lORMAL 

SHOCK  DIFFUSER 

- SUBSONI 

c 

DIFFUSER  - 

CONSTANT- APE A EJECTOR 

DATA: 

GS 

- O. 1 SSR00E+O 1 

GP 

= 

O. 1 34000E+O 1 

GM 

- 0. 1 37323E+0 1 

I-IHPMWS 

= 

0 . 1 6S375E+0 1 

A7AR 

- O. 153968E+01 

A7AS 

= 

0. 3 E 9 5 4 4 E + 0 1 

MPAS 

= 0. 577684E+01 

MR 

■ 

0. R1S000E+01 

MS 

= O.49555SE+01 

M7 

0 . 4 355R0E+O0 

PS0PR 

— 1.1.  1 ,-*8AS'~|F+H4 

PSOPRO 

= 0. 183333E+03 

T S0T R0 

= 

0. 7S137SE+O0 

P7FR 

= O. R55548E+0R 

T7T2 

0. 1S091 1E+01 

P70PRO 

= O. R74940E+O1 

T70TR0 

=: 

0. S0R0S1E+O0 
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7.4.2  CGLDOP  Sample  Case  No.  1 (Cont.) 


HOPMRL  SHOCK  DIFFUSER  - 
CONSTANT-AREA  EJECTOR  - 
DAT  Ft : 

OS  = 0. 15620OE+01 
CM  = 0. 1 37328E+01 
H7A8  = 0. 153968E+01 
R8A7  = 0.C00O00E+01 
I'lc  - 0.21 S000E  + 0 1 
NT  - O.435520E+00 
PS OP 2 = O. 193689E+04 
PS0P20  = 0. 1 83333E+03 
PS P£  = O . £804 1 5E +0£ 
P88P80  = U.278784E+81 

SUBSONIC  DIFFUSER  IiRTfl: 

PS D - 0 . 998 1 59E+00 

U-1  - 0.  1 Cm  c!c8E+0  1 

Nr  = 0. m 35520E+Q0 
PSP?  = 0. 10978 IE +01 
P80P70  = O. 9981 S9E+00 


SUBSONIC  DIFFUSER  - 
SUBSONIC  DIFFUSER 


CP  - 0.1 S4000E+01 
NNPMNS  = 0. 188375E+01 
R7R6  - 0. 369544E+01 
WPWS  = 0. 577884E+01 
MS  - O.495558E+01 
MS  - 0. 199S07E+OO 

TS0T£0  = 0. 7S137SE+00 
I oT 8 = u. 1 OjSccE+O 1 
7 80T  £0  = 0 . 8080k 1 E+00 


R8R7  --  0 . 20000OE+0 1 
MS  - 0.  199S07E-+O0 
1ST 7 = 0. 10877SE+O1 
T80T7©  = 0. 100O00E+01 


TO  RESTART  PROGRAM  ENTER  "YES" 
TO  STOP  PROGRAM  ENTER  "NO" 

YES 
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7.4 

C LG  DOF' 


POINT 
POINT 

POINT 
POINT  ~i 


.3  CLGDOP  Sample  Case  No.  2 
0 AMPLE  OUTPUT  DATA:  CASE  MO.S 


HIGH  ENERGY  CHEMICAL  LASER  SYSTEM  SIMULATION 
OHE-D I MENS  I ON AL  AHALYS I S 

A. L.  ADDY 
C.D.  MIHKELSEN 
M.P.  SANDBERG 

1 JANUAR'T  7 6 

MECHANICAL  ENGINEERING  DEPARTMENT 
UNI "EPS I TY  OF  ILLINOIS  AT  URBANA-CHAMPAIGN 
URBAN A*  ILLINOIS  61801 


SYSTEM 

LASER 

SSE  SOLUTION  FOR  MINIMUM 
DATA: 

CA'  'I TY  ENTRANCE  COND I T I ON 

P60P1 

1“ 

Ml 

- O.SS30O0E+01 

GS 

- 0.  lT'bSOOE+Ol 

PIP  10 

= O . 880 1 8SE-0 1 

riTio 

= O. a171SlE+0O 

LASER 

CA"ITY  EMIT  AND 

SUPEPSONI 

C-SUPERSON I C 

E JEC TOR  ENTRANCE  C OND I T I ONS 

GS  - U.  15bS0U£+01 

MS  - 0 . SS30O0E+0 1 ASA] 

= O. 10000OE+01 

PSP  1 

- 0. 10O00OE+01 

TST 1 

= O. 10O000E+01 

PS0P10 

= O. 10O000E+81 

T SOT  10 

= 0. 1 0000OE+0 1 

SUPERS 

ON 1C -SUPERSONIC 

EJECTOR  PRIMARY  NOZZLE 

EMIT  C 
GP 

ONDITIONS 
- 0. 1 3400OE+01 

MNP NWS 

= O.  14781 5E+01 

WPMS 

M6 

= 0 , 5897S6E+0 1 
- O. 31 SCO IE +01 

A6A1 

= 0 . i S 4 0 J 9 E + 0 1 

P60P 1 

= 0. 1S3O06E+03 

T60T 1 

- 0 . 1 8J544E +01 

P60F*  I O 

- 0. 1 0SS68E+0S 

T60T 10 

= 0 . 8073 1 3E+00 

SUPERS 

ON I C-SUPERSON I C 

EJECTOR  E 

MIT  AND 

SUBSONIC  DIFFUSER  ENTRANCE  CONDITIONS 

GM 

- U.l 0600 1 E +0 1 

MHMMWS 

= 0. 138360E+81 

WMWS 
m r 

— m . k 7c6E+0 1 

- O.47076£E+O0 

A?fi  1 

— 0. SS48J8E+01 

P7P 1 

— U . 1 057 62E+0S 

T7T 1 

- O. 1SS517E+01 

F'7uP  10 

- 0.188710E+01 

T78T10 

= 0 . 835845E+0O 
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7.4.3  CLGDOP  Sample  Case  No.  2 (Cont.) 


OLGDOP  ■-.flMPLE  IMPUT  DAT  ft:  CASE  NO.  2 

INPUT  BATA  FOP  LASER  CA'.'ITY  ANALYSIS  BY  NAMELIST. 

CURRENT  "ALLIES  ARE: 

; C fil  I 

03= 1 . 40000  j M 1 -S . O0O00  ? A2A  1 = 1. 00000 , C 0 1=0.  OOOOO i 
C 02=0 . OOOOO  j MPT S=£ 1 1 * 

4 C fti  i GS=  1 . 562  ? M 1 =5 . 0 * A2A  1 =2 . 55724  > CQ  1 = 1. 04 

INPUT  THE  "ARIAELE  TO  BE  MINIMIZED  FROM  THE  FOLLOWING  LIST 

F60P1"  DIMENSIONLESS  PRIMARY  STAGNATION  PRESSURE 
"URNS"  PR  IMAP  i ' - T 0 - S E C 0 N D A P Y MASS  FLOW  PATIO 


UF  MS 

INPUT  SYSTEM  CONSTRAINTS  BY  NAMELIST. 
CURRENT  "ALUES  ARE: 


-■  CONST  1 

F UF'  1=76.  OOUO  ? F'6uP  1 =c.'5u0 . 00  > 4 
4 C OUST 1 PSP  1=57.1 6 34  - P60P 1 =3045 . S3$ 


INPUT  FHE  EJECTOR  MODEL  FROM  THE  FOLLOWING  LIST: 


' CAE1'  CONSTANT-AREA  EJECTOR 
' SSE"  SUPERSON ' C'-SUPEF'SOM  I C EJECTOR 


C RE 


INPUT  DATA  FOR  SUPERSONIC -SUBSONIC  DIFFUSER  SECTION  BY 
NAMELIST.  CURRENT  UALUES  ARE: 

4BIFUSR 

PI  ISD=  1 . OOOOO « A4A3=  1 . OOOOO  j 4 
4 D I FIJSF:  RNSD=0 . 75  * R4h3=2 . 04- 

INPUT  DATA  FOR  EJECTOR  ANALYSIS  BY  NAMELIST. 

CURRENT  "ALUES  ARE: 

. EJECT 1 

GF  = 1 . - 0000 j MWPMWS= 1 . OOOOO  j T 60T50= 1 . OOOO0 j R8A7= 1 . O0OO0 . 
4 E JEC T 1 GF- 1 . 34 j MWPNWS= 1 . 6S375 j T60T50=0 . 76 1 376 > ASR7=2 . 


r 


i 


221 

7.4.3  CLGDOP  Sample  Case  No.  2 (Cont.) 

POINT  8 SUBSONIC  DIFFUSER  EXIT  CONDITIONS 


GM 

0. 1 3LSS1E+01 

MMMMMS 

= 

0. 1 3SJL0E+01 

IIMNS 

U . LOS 7 c'LE+0 1 

MS 

U . S 1 cj5?E+00 

ASA  1 

0. 440L70E+01 

POP  1 

o. soroosE+os 

TST  i 

0. 1O87SOE+01 

PSOP10  = 

0. 1 S70S4E+0 1 

TS0T10 

m , yj5845E+8vJ 

LASER  Cfl' 

'I TV  DATA: 

MATS 

r- 

„i 

CO  1 =■- 

0. O00OO0E+OO 

COS 

— 

0. 000O00E+00 

GS 

0.  1 T-LSOOE+Ol 

ASAl 

- 

0. 1O000OE+01 

Ml 

U. ScOOOOE+01 

MS 

= 

0. SS3000E+O1 

PSP  1 

0 . 1 OG08OE+0 1 

TST  1 

= 

O. 1O0000E+01 

PS0P10  = 

Oo 10O0O0E+01 

TS0T10 

0. 10000OE+01 

SUPEPSON I C-SUPERSON I C 

EJECTOR  DATA: 

GS  = 

0. 1 5LS00E+0 1 

GF' 

0. 134000E+OI 

GM 

0. 1 3L881E+01 

MNP M MS 

= 

0. 147S15E+01 

A7AS 

0 . SS4S  30E+0 1 

A7AL 

= 

0. 1 SO  1 03E+0 1 

MF'MS 

0 . 5 00 7 S LE + 0 1 

MS 

= 

0. SS30O0E+01 

MS 

0. 31SL01E+01 

M7 

=. 

O. 47O7LSE+0O 

POOPS  = 

0.  1 LGOOt-E+O 

PLOPS©  = 

0 . 1 0 S S L SE  -*  0 S 

TLOTS0 

= 

8. SO7313E+00 

P7P£ 

O . 1 857 LSE+OS 

T7TS 

=. 

0. 10S517E+01 

P70PS0  = 

O. 18O718E+01 

T7OTS0 

- 

0 . S 3 5 S 4 5 E + 0 O 

SUPEPSON I C - SUPEPSON I C 

EJECTOR  - 

- * 

sUBSOHIC 

DIFFUSER 

DATA: 

GS 

O. i 5LS00E+8 1 

GP 

=: 

O.  1 34000E+O1 

GM 

0. 1 3L381E+01 

MWPMWS 

— 

O. 147815E+01 

A7A£ 

O . SS4SSOE+0 1 

A7AL 

= 

0. 130\03Ef 01 

ASA7  - 

0. S00000E+01 

MF'MS 

= 

0. 5007SLE+01 

ME 

U. cSGU0UE+Ol 

ML 

= 

O . 3 1 cLO 1 E+0 1 

M7 

O. 4707LSE+90 

MS 

r: 

O. S1S55OE+00 

FLOPS  = 

O. 1S300LE+O3 

FLOPS©  = 

O. 108SLSE+0S 

TL0TSO 

- 

8. 3073 1 3E  -?-00 

PSPS 

O. SO70OSE+0S 

TSTS 

=: 

0. 10S7S0E+81 

FSOF'SO  = 

U.  ly7ot:4E+01 

TS0TS0 

= 

0 . 8 3 5 S 4 5 E + 0 O 

SUBSONIC 

DIFFUSER  DATA: 

PSD 

0. OOOO0SE+0O 

GM 

0. 1 3LS81E+01 

ASA? 

= 

0. S00000E+01 

MT 

0. 4 FO7LSE+00 

MS 

•= 

0. S1S55OE+0O 

P8P7 

0. 1 1 14S5E+01 

TST  7 

= 

O. 1 03SS7E+O 1 

P80P7O  = 

0. OO0O0SE+0O 

TS0T70 

5= 

0 . 1 00000E+0 1 

TO  RESTART  PROGRAM  ENTER  "YES" 
TO  STOP  PROGRAM  ENTER  "NO" 

NO 
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7.5  CHEMICAL  LASER  CAS  DYNAMICS  SYSTEM  PROGRAM 


7.5,1  Computer  Program  (CLGDSP) 
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7.5.1  CLGDSP  (Cont.) 
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7.5.1  CLGDSP  (Cont.) 
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7.5.1  CLGDSP  (Cont.) 
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7.5.2  CLGDSP  Sample  Input 


INPUT  HATH  FOP  LASER  CAVITY  ANALYSIS  BY  NAMELIST. 

CURRENT  UALUES  ARE: 

.fCAU 

OS-  1.400000  < Ml=  £. 000000  J He H 1 = 1.000000 

C01-  0.0000000  , COS-:  0.0000000  * NPTS=  cl  ' 

FSF‘1  = 10.00000  ! J 

4 CAM  UPTS-5'  PSPl=c0. 0 


INPUT  I HE  EJECTOR  MODEL  FROM  THE  FOLLOWING  LIST: 

"CPE"  CONSTANT-PRESSURE  EJECTOR 

"CAE"  CONSTANT-AREA  EJECTOR 

"SC E " SUPERSON I C-SURERSON I C EJECTOR 

CAE 


INPUT  THE  ITERATION  MAPI ABLE  FROM  THE  FOLLOWING  LIST: 

"MS"  PRIMARY  NOZZLE  EH  I T MACH  NUMBER 

"A7AS"  MINING  TUBE  EH I T-T O-PP I MARY  NOZZLE  EXIT  AREA  PATIO 
"WPWS"  PRIMAPY-TO-SECONDAPY  MASS  FLOW  PATIO 


WPWS 


INPUT  DATA  FOR  SUPERSONIC -SUBSONIC  DIFFUSER  SECTION  £■',  NAMELIST. 
CURRENT  i 'BLUES  APE: 


SBIFUSP 

F'NSD=  1 . 00000O 


R4h3= 


1 . 000000 


TBIFUSP 


INPUT  DATA  FOP  EJECTOR  ANALYSIS  BY  NAMELIST. 
CURRENT  UALUES  ARE: 


REJECTS 

GF‘=  1 . 400000 
NS=  1.01 0000 

ASA7  = 1 . 000000 


MHF'MWS-  i.000000  i TS0T50=  1.000000  ? 

A7AS=  10.00000  < WPWS-  1.000000 


4 EJECTS  MS-4 . 5 
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7.5.3  CLGDSP  Sample  Output 


HIGH  ENERGY  CHEMICAL  LASER  SYSTEM  SIMULATION 
G!  iE  -B I Id  iS  I ONAL  RNRLYS  I S 

h.L.  Hit  BY 
C.D.  HIH.ELSEN 
IGF.  SRI  ifiBLFG 

* Jill  lURPY  YG 

MECHANICAL  ENGINEERING  DEPARTMENT 
Ul  1 1*  CPS  I TV  OF  ILL  I HO  IS  RT  UPBANR-  CHAMPA I GH 

opbama*  Illinois  si  sen 

CRE  SOLUTION  MITH  ITEPRTIOM  ON  MF'HS 
S t'sTErl  SATA : 

POINT  1 LASER  HR' 'I TV  EUTPRNCE  CONDITIONS 

Hi  - t . EOUOOOE+Ol  GS  ~ 0 . 1 4'OOOOE -*  0 1 

P i F 1 O - 0 . 1 STSOSE+OO  T 1 T 10  - 0 . GSSS5GE+ 00 

POINT  £ LRSEP  CA'»ITY  EMIT  AND  NORMAL  SHOO  DIFFUSER 
LNTFANCE  CONDITIONS 

GS  - 0. '4C000E+01 

H£  --  O.GOOOOOC+Ol  h£R1  --  0.  100000E->  01 

F£P 1 - 0. i 00000 E +01  TET1  ~ 0. iOOOOOE+Ol 

PE'OPIO  - 0.  IGOOOOE+Oi  T £01 1 0 - O.  lOGOOOE^Cu 

POINT  3 NORMAL  SHOO  DIFFUSER  EXIT  AND  SUBSONIC 
DIFFUSER  ENTRANCE  CONDITIONS 

GS  - 0.  HOOOOE+Ol 

HI:  - 0 . 077  SSOE+OO  ASA1  - 0.  lOOOOOE+Oi 

F OF  l - 0.  m5OOO0E->  Oi  T3T1  - 0.  1SST50E-+01 

P 3 OF  1 0 - 0.  ."£0S?‘JE-+00  T SOT  10  = 0.  iOOOOOE+Ol 

POINT  4 SUBSONIC  DIFFUSER  Ei: ! I T AND  SUDDEN  ENLARGEMENT 
ENTRANCE  CONDITIONS 

GS  - 0.  1 '-OOOOE+Ol 

(•(•-  --  0. SrrSj0E  + 00  H4A1  - LG  1OOO0OE-+01 

F'jpi  - ft.  4c.'O0OOE->-01  T-T 1 - 0.  lSS750E-f  OI 

R40P1O  - O.7S0S74E+0O  T40T10  - 0. IOOOOOE+Ol 
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7.5.3  CLGDSP  Sample  Output  (Cont.) 


POINT  5 CONSTANT-AREA  EJECTOR  SECONDARY  NOZZLE  EXIT 
CONDITIONS 

GS  = 0. 14O0O0E+O1 

M5  - 0. O46047E+00  A5A1  = O . SSbCSSE+OO 

P5P1  ~ 0 . 4c603bE+0 1 T5T 1 ~ 0.  1 ‘*■6 1 3cE+U 1 

PC OP 10  = O.720874E+00  TC0T10  = @. lOb^OOE+Ol 

POINT  b CONSTANT- AREF  EJECTOR  PRIMARY  NOZZLE  EXIT 
CONDITIONS 

GP  = O.H0000E+01  MUF'MMS  = 0.  1O0O0OE+01 

WPWS  - 0.£65889E+01 

Mb  - 0. 450OO0E+01  AbAl  = O. 104O65E+0O 

PS0P1  = O. 19b£l£E+04  TbOTl  = O. 1SOO00E+81 

POOR 10  = O. 2507S8E+03  TOOT 10  = O. 100O0OE+01 

POINT  7 CONSTANT- AREA  EJECTOR  EXIT  AMD  SUBSONIC 
DIFFUSER  ENTRANCE  CONDITIONS 

GN  - 0. 140000E+01  N NUNNS  = O. lOOOOOE+Ol 

WMWS  - 0. S65889E+01 

M7  - O.4b£037E+00  A7A1  = O. 1O4O85E+01 

R7P1  = 0. 199966E+02  T7T1  = O. 17£6£9E+01 

F7OP10  = 0. ZS583SE+01  T70T10  = 0. 1Q6000E+01 

POINT  S SUBSONIC  DIFFUSER  EXIT  CONDITIONS 

GN  - O. 140808E+01  HUMMUS  = O. lOOOOOE+Ol 

UMUS  - 0. 385889E+01 

MS  - 0.462037E+O0  ASA 1 - O. 1O4O85E+01 

F-8P1  0.  1999b6E+0£  TST1  = O.  17£b£9E+01 

P80P10  = O . 2S5S3SE+0 1 T80T1O  - 0. 1O0000E+01 

LASER  CAN I TV  DATA: 

NF'TS  - 5 

CO.l  - O.000000E+G0  COS  - 0 . O000OGE  + 00 

GS  = 0.  14O00OE+01  AS'A  1 - 0.  1OO000E+01 

Ml  - 0.  S00000E+01  MS'  - O.800000E+01 

F'SF'1  = 0.  1O0000E+01  TST  1 - 0.  100000E+01 

PS0P10  = 0.  100000E+01  TS'OT  10  = 0.  100000E+01 

NORMAL  SHOCK  DIFFUSER  DATA: 

GS  = 0. 14000OEf01  RNSD  = O . 18000OE+81 

MS'  - 0.  S00000E+01  M3  - 0. 57735OE+00 

PSPS  = 0.45000OE+01  T3TS  - 0. 16S750E+01 

P30PS0  = 0. 7S0S74E+00  T30TS0  = 0. 1800O0E+01 
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7.5.3 

POINT 


CLGDSP  Sample  Output  (Cont.) 


SUBSONIC:  DIFFUSER  EXIT 


CONDITIONS 


(Continued) 


SUBSONIC  DIFFUSER  DATA: 


PSIl  - 0.  100000E+01 
OS  = O. 14O000E+01 
M3  = O. 57735OE+00 
P4P3  =.  0.  10000OE+01 
P40P30  = 0. 100000E+01 


R4R3  = 0.10O00OE+01 
M4  - 0 . j77  350E+00 
J4T3  - 0 . 1 0 0 0 0 0 E + 0 1 
T40T3O  = 0. 100000E+O1 


CONSTftNT -AREA  EJECTOR  DATA: 


os 

~ 0. 3 4O000E+01 

OP 

ON 

- 0. 14000OE+01 

MWPMWS 

= 

A7H6 

= 0.  1 00OO0E-+  02 

WPWS 

HO 

- 0 . 45O0O0E+0 1 

M7 

— 

P OuF‘50 

= O.347S6OE+03 

TO0T50 

P7P50 

— 0 . OS4522E+0 1 

T7T50 

— 

P70P50 

--  0.4103SSE+01 

T 70 I SO 

NORMAL 

SHOQ  DIFFUSER 

- SUBSONIC 

CONST ANT -AREA  EJECTOR  DATA: 


OS 

= 

o. 140000E+01 

OP 

1 

ON 

= 

0. 140000E+01 

MWPMWS 

55 

A7A2 

0. 1 04OO5E +0 1 

A7A0 

‘JT 

HP  NS 

- 

0. 2O5SSOE+0 1 

M2 

MO 

- 

0. 4500O0E+01 

M7 

PO0P2 

0.  l:‘O212E+04 

P00P20 

- 

O . 2 j07 OSE+03 

TO0T20 

P7P2 

0.  r?“?'?OOE+02 

T7T2 

P7oP20 

- 

G. 295S38E+01 

T70T2O 

- 

NORMAL  SHOD  DIFFUSER  - SUBSONIC 
CONSTANT-AREA  EJECTOR  - SUBSONIC 


DATA! 

OS- 

0.  14O0O0E+01 

OP 

ON 

- 

0. 140000E+01 

MWPMWS 

A7A2 

0. 1O40O5E+01 

A7A6 

ASA7 

z. 

0. 100000E+01 

WPWS 

M2 

zr. 

0 . 200000E+0 1 

MO 

M7 

z. 

0. 4O2037E+80 

MS 

PO0P2 

- 

0. 13O212E+04 

PO0P20 

= 

0 . 25 0 7 0 S E + 0 3 

TO0T20 

PSP2 

0. 1339OOE+02 

T ST  2 

PS0P20 

= 

0. 235838E+01 

T 80T20 

0. 140000E+01 
0. 100000E+01 
0.  £C?SSSE-r01 
0.  •*S£037E+00 
0. 10O000E+01 
O. 353053E+00 
O. 1OO000E+01 

DIFFUSER  - 


0.  H0000E+01 
0. 100000E+01 
0. 100000E+0C 
0 . 30O000E-+  0 1 
0. 4O2037E+00 

0. 10OO00E+01 
0.  172O29E+01 
0. 100O00E+01 


DIFFUSER  - 
DIFFUSER 


0. 140000E+01 
0. 1000O0E+01 
0. 1Q0O00E+0C 
0. 265889E+01 
0.4  5 0 0 0 0 E + 0 1 
0 . '-rSc037E  + 00 

0. 10000OE+01 
0.  172029E+01 
0. 100000E+O1 


279 


7.5.3  CLGDSP  Sample  Output  (Cont.) 

POINT  S SUBSONIC  DIFFUSER  EXIT  CONDITIONS 


SUBSONIC  DIFFUSER  DATA: 

PSD 

= 0. 100000E+01 

GM 

= 0. 140090E+01 

AoA7 

N7 

= O.46SOS7E+00 

MS 

PSP? 

= 0. 100000E+01 

TST7 

PS0P7O 

= 0, 100000E+01 

T o0T  7 

:EST ART  PROGRAM 

ENTER  "YES" 

TO  STOP  PROGRAM  ENTER  "NO" 


(Continued) 


= 0. 100000E+01 
= O.  46SO37E-+00 
- O. 100000E+01 
0 = O. 100000E+01 


END  OF  EXECUTION 
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7.6  AN  EXPERIMENTAL  INVESTIGATION  OF  THE  PERFORMANCE  OF  A CONSTANT- 
AREA,  SUPERSONIC-SUPERSONIC  EJECTOR  FOR  VARIATIONS  IN  THE  MIXING 
TUBE  LENGTH- TO-D I AMETER  (L/D)  RATIO 

The  results  presented  in  Section  3.0  of  this  report  were  prin- 
cipally for  variations  in  primary  and  secondary  stream  Mach  numbers, 
Mpj  and  Mgi  , and  secondary-to-primary  area  ratio,  The  ef- 

fects of  variations  in  the  mixing  tube  length-to-diameter  ratio  were 
investigated  only  indirectly  by  means  of  the  measurement  of  the  flow 
non-uniformity  at  the  mixing  tube  exit  and  the  recompression  pres- 
sure rise  along  the  mixing  tube  wall.  The  need  to  investigate  the 
effects  of  length-to-diameter  ratio  was  recognized,  and  a series  of 
experiments  were  conducted  for  an  ejector  configuration  investigated 
in  Section  3.0.  This  configuration  was  defined  by:  Mgi  = 2.00,  Mp x = 
2.50,  and  A^/A^  ~ 0.88;  experimental  data  were  presented  for  this 
ejector  in  Figs.  3.3  - 3,7,  and  9.  For  this  series  of  experiments, 
the  experimental  procedure  was  the  same  as  outlined  in  Section  3.2, 
and  the  mixing  tube  ler.gth-to-diameter-ratios  were:  L/D  = 5,  7.5,  10, 

12.5,  and  15. 

The  results  of  this  investigation  are  presented  in  Figs.  7.6  - 
1,2,  As  is  evident  from  these  figures,  the  overall  ejector  perform- 
ance was  nearly  identical  for  all  ejector  configurations  with  L/D>^ 

7.5.  The  largest  values  of  Pf0/Psl  versus  Wp/Wg  were  found  for  a mix- 
ing tube  L/D  = 10  (the  configuration  investigated  in  Section  3.0). 

The  values  of  PfiC/Psl  versus  Wp/Wg  for  mixing  tube  L/D  =7.5,  12.5  and 
15  were  slightly  less  than  those  for  a mixing  tube  L/D  = 10.  Ejector 
performance  at  a mixing  tube  L/D  = 5 was  significantly  degraded  rela- 
tive to  the  mixing  tubes  with  large  L/D  ratios.  This  poor  perform- 
ance is  due  to  the  inadequate  mixing  and  recompression  length  within 
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the  L/D  = 5 mixing  tube.  For  the  ionger  mixing  tubes,  frictional 
losses  may  tend  to  reduce  the  performance  more  than  the  gains  made 
in  achieving  a more  uniform  velocity  profile  at  the  mixing  tube  exit; 
compare,  for  example,  the  performance  results  for  L/D  = 7.5,  10, 

12.5,  and  15  in  light  of  the  known  flow  non-uniformities  for  an  L/D  =>  10, 
Figs.  3.3  - 7,9. 


Compres 

Diameter 
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